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SUMMARY AND CONCLUSIONS 


1. The very high lightning flash density in many parts of the United 
States and the large size of the SPS rectenna require us to incorporate 
lightning protection systems in the rectenna design. 

2. A distributed lightning protection system is described in this report 
that will protect the rectenna components from direct lightning strike 
damage and will, in addition, provide reduced induced lightning effects 
in the power and control circuits. 

3. The proposed lightning protection system should be incorporated as a 
structural member of the rectenna support system; viewed as such, the 
lightning protection system will not appreciably increase the total 
material requirements for the rectenna unless materials are used that are 
incapable of safely conducting lightning currents. 

4. The lightning protection design places the conducting elements so that 
the microwave shadow cast by protection systems falls along the upper 
edge of the billboard on which it is mounted (and the lower edge of the 
next billboard to the north); these shadow areas are only a slight 
fraction of the collecting area, so the protection elements produce very 
little, if any, additional power loss to the rectenna as a whole. 

5. Individually the microwave diodes are self-protecting with respect to 
"average" lightning and those near the center of the rectenna are safe 
from extreme lightning. However, the series connection of the diodes to 
form 40,000 V strings creates a protection requirement for the string. 
Standard surge protection practices are necessary for the string. 

6. Electric power industries usually attribute 10% of the cost of power 
transmission equipment to lightning protection requirements. If this 
factor is not already included in cost estimates, it should be added. 
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SUMMARY OF THE RECOMMENDED LIGHTNING PROTECTION DESIGN 


Based upon our research, computer simulations, and laboratory tests with 
a scale model, we recommend a distributed lightning protection system that 
employs a horizontal conducting member with points and grounds placed at every 
bay or billboard (14.69 meters apart). This conf i guration not only provides 
greater protection than other configurations that were evaluated, it is more 
easily integrated into the structural design of the rectenna. The recommended 
system is shown in Figure 1. 
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PREFACE 


The objectives of this study are to evaluate the hazard posed by light- 
ning flashes to ground on the SPS rectenna and to make recommendations for a 
lightning protection system that will provide sufficient protection to the 
rectenna. For purposes of this study, the SPS rectenna design is based upon 
the data supplied to us by Rockwell International in July, 1978. 

This study has four major components, each with several elements of 
investigation. The components were: lightning distribution; lightning 
interactions; rectenna damage estimates; rectenna protection. The elements of 
each component are listed in Table A. The study plan was to proceed from top 
to bottom evaluating the elements listed in each component; work proceeded in 
a parallel manner for the four components. The organization of this final 
report reverses this order by presenting the more important results of the 
study first, then following this with the material and considerations leading 
to the conclusions. 


vii 



TABLE A 


Rectenna Electrostatic Protection 


1 . 


2 . 


3. 


4. 
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I. ANALYSIS OF LIGHTNING ROD PROTECTION CAPABILITIES FOR A CONFIGURATION 
SUITABLE TO THE RECTENNA 

1. Cone of Protection Considerations: 

I. 1.1 Definition and Considerations 

The capability of a vertical conductor to attract a lightning flash 
is described by the cone-of-protectio n, or perhaps more accurately the 
cone-of-attraction. In theory, any lightning flash that would have 
entered this cone had the vertical conductor not been in place, will 
strike instead the conductor and be shunted to the ground. The method by 
which this process takes place is as follows: 


The lightning stepped leader creates high voltages over a wide 
area on the rectenna because of the large charge on the leader 
tip. At points on the rectenna where the electric field reaches 
breakdown values due to local enhancement factors, upward 
propagating sparks are initiated which move to meet the downward 
propagating stepped leader. The upward propagating spark which 
first makes contact with the leader completes the electrical circuit 
and the lightning flash current will pass through the structure that 
initiated the successful upward going spark. 


The cone of protection is primarily a function of the height of 
the vertical conductor because of the field-enhancement factor which 
enables the taller object to initiate the upward spark before lower 
objects. Other factors enter into the consideration of the cone of 
protection, such as the charge on the leader tip and the velocity of the 
leader, because these factors strongly influence the timing of the 
production of upward sparks and the height at which the spark and leader 
meet. In general, the results of research into this subject have shown 
that the larger the leader charge, then the larger the angle 3 of the 
associated cone of protection. Since larger leader charges are usually 
associated with the larger lightning currents, we find a fortunate result 
that the cone of protection increases with the potential hazard of the 
lightning flash. 

It follows then that the angle 3 of the cone of protection (See 
Figure 2) varies with the particular lightning flash. 3 = 45° is a very 
commonly used design angle in the United States and many of the examples 
in this report employ 3= 45°. 
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Figure 2 


1.2 Distributed Lightning Protection Systems 

The cone of protection and the experimental data used to evaluate 
are specifically related to the single elevated point, and in most cases 
the system under consideration is 10 to 100 meters in height. As will be 
seen later, lightning protection of the rectenna falls into a class of 
structures that requires distributed lightning protection tactics. 

Figure 3 illustrates a distributed system used by power transmission 
companies. The main point is that the cone of protection concept is of 
limited usefulness in the total protection problem. We will use it on 
the panel and billboard scale as a technique to make a comparative 
assessment of capabilities of various configurations. 

2. Lightning Rod Protection Configurations Compatible with the SPS Rectenna 

We have considered three different configurations of lightning rod 
systems in this effort. In the smallest scale system considered each rectenna 
panel (0.74m in width) had a short lightning rod attached; see upper example 
in Figure 4. In the medium scale system each rectenna support structure 
(14.69m apart) or billboard will have an attached lightning rod; see middle 
example in Figure 4. And, in the distributed protection system, short termi- 
nals located on each rectenna support structure (14.69m apart) were connected 
by horizontal conducting structures; see lower example in Figure 4. 

As seen in the analysis of the billboard scale system, it is impractical 
to seriously consider larger scale systems. 
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POWER LINES EMPLOY DISTRIBUTED LIGHTNING PROTECTION SYSTEMS. THIS- ILLUSTRATION 
SHOWS A "STATIC" OR GROUNDED PROTECTION WIRE TAKING A STRIKE AND PROTECTING THE 
POWER LINES BELOW. 


FIGURE 3 
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PANEL SCALE LIGHTNING PROTECTION SYSTEM 

FIGURE 5 



2.1 Lightning Rod Protection at the Panel Scale 

In this system configuration a relatively short lightning rod is 
positioned at the top of each panel and oriented perpendicular to the panel 
face (see Figure 5). Conceptually the rod is centered on the top of the 
panel, but in practice it could be on the panel edge without altering the 
results of this analysis. 

Let a be the inclination of the rectenna. Figure 6 illustrates the case 
where 0, the angle of the cone of protection, is greater than a. This figure 
applies only to the conditions in the vertical plane that passes through the 
lightning rod and is perpendicular to the rectenna face. In this particular 
projection it appears that the short (example 0.74m) lightning rod on the 
panel provides adequate protection to the rectenna. In other projections we 
see that there are, however, "holes in the armor. 

Figure 7 is an enlargement (xlO) of the lightning rod portion of Figure 
6, and defines the parameters to be used in the following discussions. The 
cone of protection intersects the plane of the rectenna to form conic 
sections : 

(1) If o+ 0 = 90 0 the intersection is a parabola. 

(2) If a + 0< 90° the intersection is an ellipse. 

(this is the case illustrated in Figures 6 & 7) 

(3) If a+ 0 > 90° the intersection is a hyperbola. 


If we now look at the intersection of the cone of protection with 
the panel for the specific cases above, we see the emergence of the 
protection problem with this type of lightning rod protection configura- 
tion. From the geometry of Figure 7 we see that the axis of the cone is 
at l = L tan a and that the vertex of the conic is at d = L 
tan (0- a) relative to the top of the panel. 
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In a coordinate system defined in the rectenna plane with the origin 
at the axis of the cone and the y axis directed north (toward top of 
rectenna) and the x axis directed east, the equation for conic is: 


2 2 2 2 2 2 2 
x cos a y (cos a - sin a tan 8)cos'a 2y sin 

2 2 2 2 
L tan 6 L tan B 

For the parabolic solution this equation reduces to: 


a cos 
L 


a 


1 


2 

x 


2L sin 2 g , L 

cose COS a ^ 2 COS 6 COS a 


In figure 8 we have plotted the intersection of cones of protection for 
three lightning rods of lengths 0.185m (= V 4 panel width), 0.37m (= V 2 panel 
width), and 0.74m (= panel width.) 


In these examples the rectenna inclination angle a is taken to be 45° and 
the cone of protection 8 is equal to 45°. The resulting intersections are 
parabolas for the cases depicted in Figure 8. For the parabolic solution the 
cone of protection is parallel to the face of the rectenna in the vertical 
plane bisecting the panel (The view of Figure 6 and 7 except that 
here a = 8 = 45 u ) . 

At lower latitude sites (below 40°) the rectenna inclination angle a is 
less than 45° and the 45° cone of protection intersection becomes an ellipse; 
in Figure 6 the vertical projection illustrates the intersection in the plane 
through the lightning rod. The elliptic solutions leave regions along the 
base of the rectenna unprotected. Hence, the parabolic solutions of Figure 8 
and the table (Fig. 9) represent maximum protection capabilities of the cone 
of protection with the panel scale protection conf iguration. The small 
ellipse in Figure 11 shows the cone of protection intersection 
for a = 40 °, 8 = 45 °, and L = 0.74m. 

2.2 Lightning Rod Protection at the Bay or Billboard Scale 

In this system a longer lightning rod is placed at the center (or 
end) of each bay or billboard making them 14.69m apart. The mathematical 
description here is identical to that for the panel scale system (2.1). 
Only sizes are different. Figure 10 illustrates the billboard scale 
system. 
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FIGURE 8 


THE INTERSECTION OF THE CONE OF PROTECTION WITH A RECTENNA 
PANEL (THE CURVED LINES) SHOWN IN THE PLANE OF THE PANEL, 
LIGHTNING ROD LENGTHS = kj k AND 1 TIMES THE PANEL WIDTH 
ARE SHOWN PROJECTED VERTICALLY ONTO THE PANEL, 



PARABOLIC TYPE SOI UTIONS 


UNPROTECTED AREA 

ROD LENGTH IN MFTFRS UNPROTECTED AREA IN % X ENHANCEMENT FACTOR 


,185 1.1% 2,9% 
.37 ,55% 1,5% 
,7H .28% .7HI 


FIGURE 9 





To illustrate the cone of protection concept for this configuration we 
use as an an example, a = 40 °, 0 = 45 °, and L = 7.35m ( = 1/2 billboard 
width). The resulting intersection is a portion of an ellipse and is shown on 
Figure 12. Even if these long (7.35m) lightning rods were placed every 
14.69m, a significant fraction of the rectenna (6.7% or when weighted by 
enhancement factor 18%) is unprotected (i.e. is not inside a cone of protec- 
tion) . 


Furthermore, there are serious mechanical problems associated with sup- 
porting these long (i.e., over 22 feet) lightning rods. We think these 
examples are sufficient to demonstrate that configurations employing fewer 
lightning rods at longer spacing decreases protection and creates structural 
problems that ultimately will increase the total materials requirement. 

For example, if we were to increase the length of the lightning rod in 
this configuration to the point that it could offer protection to the 
billboard in front of the one on which it is mounted ( i.e. to the south), 
then with the appropriate phasing of rods between rows of billboards we could 
get total protection in the cone of protection context. The length of the 
rods would need to be 12m in order to provide this coverage. 

2.3 The Distributed Lightning Protection System 

The distributed lightning protection approach replaces the many lightning 
rods with a continuous horizontal conducting structure, as depicted in Figure 
13. The region of protection now becomes the volume beneath two planes whose 
intersection is the horizontal protecting structure. This protection tactic 
is essentially the one employed by the power transmission companies. The 
angle between the protecting planes and vertical is variable; 45° is thought 
to be adequate but some designs use 30° for extra protection. This line is 
called the "static" by the power companies and this term is used here for 
convenience. 

Figures 7 and 8 provide the correct geometric considerations for the 
distributed lightning protection if we interpret the end point of the 
lightning rod to be the location of the static. We note that the figures 
apply anywhere along the rectenna, not just in the specific locations required 
by the lightning rod analysis. 

For consistent comparisons with the other lightning rod systems we will 
use a = 45°. Since a < 45 0 for rectennas below 40° latitude, the top edge of 
the rectenna is protected by the static for any value of L, the displacement 
distance. If we try to use a smaller, more conservative value for 0, we will 
run into problems in protecting the top edge of the rectenna with any system 
tht does not cast a radio shadow on an active rectenna surface. The design 
constraint that we will use to specify L will be that the southward plane of 
protection intersect the rectenna surface at the base. Therefore, 

L = 12.2m tan (45 0 - a) . 

For a in the range 45° to 30°, L has the range of values 0m to 3.3m. This 
simple analysis ignores the protecting capability of the immediate southward 
row of the rectenna on the base of the row being considered. When these 
additional protective effects are considered we find that: 
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L = 6.1m (1 - tan a) 

For a in the range 45° to 30°, L now has the range Om to 2.6m. 

Figure 13 gives the configuration of the distributed lightning protection 
system for a = 30°, which represents the most difficult situation to 
protect. In this situation the static is displaced by 2.6 meters from the top 
edge of the rectenna; note that the 45° planes of protection provide total 
coverage of the rectenna. 

We wish to emphasize that the set of horizontal statics not only provide 
total protection in the sense that lightning flashes are expected to hit the 
statics instead of the active rectenna surfaces but that this system also 
reduces the induced voltages and currents in the rectenna. We estimate that 
induced charges, currents, and potentials are reduced by 1/2 by the static 
protection system. 
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FIGURE 12 

DISTRIBUTED LIGHTNING PROTECTION SYSTEM 




FIGURE 13 


DISTRIBUTED LIGHTNING PROTECTION SYSTEM 
ILLUSTRATING FORWARD AND BACKWARD PRO- 
TECTION FOR SMALL INCLINATION ANGLES 



II. SIMULATIONS OF LIGHTNING STRIKES TO THE SPS RECTENNA WITH AND WITHOUT 
PROTECTION 


A series of experiments were performed in our electrostatic test chamber 
with a scale model of the SPS rectenna. The experiments consisted of exposing 
the model rectenna to a series of high voltage discharges produced with a 
Tesla coil . 

The strikes to the rectenna were photographed using time exposures in a 
darkened room. A wire from the upper plate conducted the discharge to the 
vicinity of the model rectenna and provided us with a limited control over the 
area of the strike. This allowed us to keep the strikes near the volume in 
focus by the camera. 

Different areas of the model rectenna were protected by different 
systems, and one area was unprotected. The following paragraphs describe 
samples of these experiments: 

1. The Unprotected Rectenna 

Most of the strikes were to the upper edge of the billboard because of 
the larger enhancement factor at that point. Several strikes to the billboard 
face occurred. 

In Figure 14, we see two strikes to the unprotected billboard section, 
one of which is to the billboard face. Notice that these strikes are perpend- 
icular to the face when near the face; we would anticipate this because the 
equi potential lines are nearly parallel to the face here. 

In Figure 14, we also see for comparison the' three lightning protection 
systems modeled. To the left is the billboard scale system; to the right is 
the panel scale system; and behind the flashes is the distributed lightning 
protection system. 

2. The Panel -Scale Protection System 

The next three figures are examples of strikes photographed on the 
section of the model rectenna that was protected by the panel -scale lightning 
protection system. 

In Figure 15, we see two strikes on the same billboard, both of which 
terminate on the panel-scale lightning rods. 

Figure 16 shows two strikes from a different view going to two different 
billboards. The panel-scale protection system here is seen to protect only 
the front billboard. Protection is probably greater for real lightning 
because in our experiments we artificially bring the "leader tip" very close 
to the billboard with the wire. 

Multiple strikes to the panel-scale protection system are seen in Figure 
18. One of the strikes goes directly to the billboard face, this type of 
failure will occur in nature, but with lower probability than illustrated 
here. 
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3. The Billboard-Scale Lightning Protection System. 

Two sets of experiments were made with the bill board- scale lightning 
protection system. The one illustrated in Figure 19 corresponds to rods of 
length 7.35m. (A second series of strikes were made with rods cut to one- half 
of this length, but these were photographed in color and are not suitable for 
this report.) Figure 19 illustrates the capability of these long rods to 
direct lightning to the desired point. 

In Figure 20, we have a side view of a bil 1 board- seal e protector taking a 
strike and protecting the billboard- face. Figure 21 illustrates the "hole in 
the armor" of the bil 1 board- seal e lightning protection system. Two flashes 
strike the protection system, but a third strikes the billboards between two 
protectors, as predicted in Figure 12. With real lightning this is less 
likely to happen, but it can and will occur. 

4. The Distributed Lightning Protection System. 

The displacement distance of the static from the billboard was scaled 
from 0.74m to make it correspond to the height of the panel -scale protection 
system. Fewer fail ures-to- protect were observed with this system but they did 
occur. With real lightning, they would be even less likely to occur. 

In Figure 22, we see two strikes to two different billboards from the 
side view. Figure 23 shows two strikes to the same billboards, which were 
rovided with a distributed lightning protection system. One strike is to the 
terminal support rod at the billboard edge, which is the preferred point of 
strike. The other strike goes to the horizontal static line between the 
terminal support rods. 
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FIGURE 14 






FIGURE 16 



FIGURE 17 





FIGURE 19 




FIGURE 21 



III. GROUNDING CONSIDERATIONS FOR THE PROPOSED LIGHTNING PROTECTION SYSTEM 


The thundercloud charges induce a large surface charge on the rectenna 
below the cloud; during the stepped leader period even larger surface charges 
are induced on the region below the leader tip. Most of the current flowing 
during the return strokes of the lightning flash must be distributed by the 
grounding system to connect with the induced surface charges. If adequate 
paths for these currents are not planned and provided, the lightning will make 
its own paths. Most of the induced surface charge will reside on the horizon- 
tal statics of the recommended distributed lightning protection system. The 
primary grounding system described here is to provide low impedance paths for 
the redistribution of the induced surface charges and the part of the light- 
ning charge that resides on the rectenna surface. 

1. Primary East-West Grounding 

It is absolutely necessary that the horizontal statics have a good low 
impedance connection at billboard edges. The static should appear to be a 
continuous very low impedance conductor in the east-west direction, as illus- 
trated in Figure 24. 

2. Primary North-South Grounding 

It is also necessary that the statics are mutually grounded in the north- 
south directions; there are two methods of achieving this: 

2.1 Periodic connections north-south at the level of the statics. If 
these north-south statics are aligned along the billboard edges, then 
there will be little power loss due to microwave shadows (See Figure 24.) 

2.2 Interconnect grounding in the north- south direction at the surface 
or sub-surface level (see figure 25) can also be used, but this approach 
creates a higher impedance to north-south currents on the static system. 

2.3 A surface level grounding network is required in addition to the 
primary static grounding network. The surface network must handle the 
redistribution of induced charges on the rectenna surfaces and power 
distribution systems and it provides a safe working environment at the 
surface level. East-west continuity with low impedance connections must 
be provided at the base of the rectenna support structures, and north- 
south continuity with low impedance connections as discussed in 2.2 and 
illustrated in Figure 25 must be provided. Figure 26 highlights the 
surface level grounding network. 

2.4. Interconnections between the primary and surface grounding networks 
should be provided by the vertical conductors located at every billboard 
upper corner; these are the same structures on which are mounted the 
terminals and supports for the statics. The vertical interconnections 
are highlighted in Figure 27. 

2.5 The ultimate or final component of the grounding system is the tie- 
in to Earth ground. At regular intervals in the rectenna a deep earth 
qroundi nq rod must be driven into, .the soil to make good contact with a 
conducting soil for earth grounds 
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The organization of the earth grounding system should be along diagonals, 
as illustrated in Figure 28. Here we see that the placement of earth ground 
at every fourth billboard but on a diagonal produces a grid such that 
lightning striking the primary grounding network will never have to travel 
more than 30 meters along the east-west conductors before finding a ground, or 
32 meters along the north- south conductors (for a rectenna with a 40° inclina- 
tion angle). 
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RECTENNA LIGHTNING ROD SYSTEM 



FIGURE 25 
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IV. MATERIALS AND SPECIFICATIONS FOR LIGHTNING PROTECTION 


It is premature to specify the final form for the materials for the 
lightning protection system. We think that the system should be integrated 
into the structural design of the rectenna itself; in this case many other 
considerations are necessary in addition to the capability to conduct 
lightning currents. The data displayed in Figure 29 (H. Baatz, Protection of 
Structures, in Lightning Vol . 2, ed. by R.H. Golde) is useful for order-of- 
magnitude estimates of the lightning current requirements. 

Example: If the design permits a 100° C temperature rise in an aluminum 

member carrying lCr Amps for lCr seconds, we need approximately 3 mrrr crossec 
tional area of aluminum material in the conductor. ~Note that the recommended 
crossections for building codes are larger 80 m) ) indicating designs for 
lower temperature operation plus safety margins. 

The lightning conductor need not be solid. From a structural point of 
view a tubular or other extruded shape would be preferable. Such configura- 
tions are compatible al so with the lightning protection recommendations. 
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Specific values of materials for wire 


Material 

Steel 

Copper 

Aluminium 

Density (g/cm -3 ) 

7-7 

8-92 

2-7 

Electrical resistance (ft mr 1 m _1 ) 

0-17 

0-0178 

0-029 

Heat (cal “C" 1 g- 1 ) 

0-115 

0-093 

0-023 

Melting point (°C) 

1,350 

1,083 

658 



Temperature rise of conductors as function of current square impulse per 
cross-section square; Cu = copper, A1 = aluminium, ST = steel, CRS = cor- 
rosion-resistant steel. 


Cross-section for lightning conductors 





Dimension 

Installation 


Cross-section 

Rod 

Strip 

components 

Material 

(mm 2 ) 

(mm, radius) (mm x mm) 

Air termination 

Steel, galvanized 

50 (25)“ 

8 

20x2-5 

Rods up to 0-5 m 

Steel, stainless 

110 

12 

30x3-5 

long 

Down conductors 

Copper 

50 (16)“ 

8 

20x2-5 

Conductors in 

Aluminium*’ 

80 (25)“ 

10 

20x4 

ground 

'"Steel, galvanized 



0-5 mm 

Sheet metal 

Copper 

Aluminium, Zinc 



0-3 mm 
0-7 mm 


[Lead 



2.0 mm 


“ Lowest cross-sections used in some countries. 
b Not for use below ground. 


FIGURE 29 
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V. ESTIMATE OF POWER LOSS FROM THE BEAM 


A rough maximum estimate of the power loss from the microwave beam due to 
the lightning protection devices can be obtained by assuming that the micro- 
wave shadow cast by the static lightning protection system is twice the 
crossectional area of the devices. We assume that the conductors are 2 cm 
wide of 1 mm thickness tubular material, providing 60 mrrr of crossectional 
area for conducting. The assumed shadow of these structures is approximately 
0.6% of the rectenna area (see Figure 30.). This is a maximum estimate of the 
loss. 
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VI. MICROWAVE DIODE FAILURES DUE TO INDUCED CURRENT TRANSIENTS 


The 25 W S GaAs diodes used in the design of the SPS rectenna have not 
been produced and no failure data is available for these devices. In order to 
obtain estimates of failure power of the diodes in the design, we used the 
specification data for the HP5082-2824 microwave diode and scaled the charac- 
teristics to 25 W using the ."Wunsch relationship" described in the references 
below. We also obtained advice directly from Dr. D.C. Wunsch regarding the 
extrapolated power failure current. 

1. Defense Department Report D224-13042-1 EMP, Susceptibility of 

Semiconductor Components, dated September, 1974. 

2. Defense Department Report D224-10022-1 EMP, Electronic Analysis 

Handbook, dated May, 1973. 

3. Defense Department Report D224-10019-1 EMP, Electronic Design 

Handbook, dated April, 1973. 

Figure 31 shows the predicted failure power for 25 watt diodes, as a function 
of pulse width. 
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PULSE WIDTH - MICROSECONDS 


FIGURE 31 
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VII. COMPUTER SIMULATION OF ELECTROSTATIC FIELD AROUND AN SPS RECTENNA 


The electrostatic fields produced by the charges on the lightning channel 
induce charges on the rectenna and on the lightning protection conductors. 
Changes in this electrostatic field require a redistribution of charge on the 
rectenna system; the resulting currents can cause diode failure even with a 
lightning grounding system in place. One output of the computer simulation of 
the electrostatic field around the SPS rectenna is an evaluation of the 
induced current on the rectenna with and without the recommended lightning 
protection equipment. 

An additional output from the computer simulation is the potential around 
the rectenna billboard enabling us to estimate the enhancement factors of the 
electric field due to the billboard shape. 


The algorithm used in the simulation computes an array of values for the 
potential around the middle of five infinitely long billboards. We assume 
here that the contribution to the local potential from billboards further away 
is ignorably small. The surface charge distribution on the billboards is 
simulated with ten infinitely long line charges evenly spaced along the bill- 
board. The value for the line charges is determined interact! vely with the 
computer to produce a zero potential contour that has the same shape as the 
billboard. Figure 32 illustrates this simulation. 


In order to compute the potential, we will need U(x,y), the electrostatic 
potential at a point (x,y) in free space, where the coordinate system is such 
that the line of electrical charges giving rise to the potential is located at 
the origin. If we call the y-coordinate the height h, then U(x,H) is the 

at x and h of a line charge A (coulomb/meter) at a 
x = 0. There is also a contribution to U 


electrostatic potential 
height d directly above the point 
from the image charge. Thus, 


U(x,h) = - 


2 ire 


In £ 


+ (h - d)‘ 


1/2 
] • 


x“ + (H + d) “ 
around the rectenna 


may be calculated. 


From this, the potential distribution 
Let U(l,h) be the potential at x = 1 and y = h due to a periodic system of 
line charges simulating the rectenna (see Figure 31.) We then have that 

w m n 0 - L[i - 1] - X.f + (h - sX. f 1/9 

U(l,h)= i= \ - - ^- o r ^U,] 172 , 

11 J 1 ^ o (1 - L[i - 1] - X. f + (h + sX, f 

J J 

where the free-space value for the dielectric constant is assumed and where 


i = Billboard number, 
j = Line charge number on billboard i, 
s = Slope of billboard (= tan a), 

M = Number of line charges (= 10), 

N = Number of billboards (= 5). 
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In the presence of a uniform electric field of 100,000 volts/meter 
(directed upward), ten line charges have been selected to produce the array of 
values shown in Figure 33. Three potential contours have been sketched (zero, 
10,000 V, and 100,000 V) around the ten line charges on the billboard. The 
zero contour follows closely the position of the billboard surface, as 
required by the simulation algorithm. Note how closely spaced the contours 
are at the top edge of the billboard. Electric field enhancement factors of 
at least 6.5 exist in this region based upon our simulations. Higher 
resolution simulations would be required to refine the enhancement factor 
estimates. 

The values obtained for the 10 individual line charges found for the 
solution shown in Figure 33 are (in pCoul./m): 

0.36, 0.465, 0.572, 0.679, 0.924, 1.02, 1.14, 1.78, 2.91, 4.14. 

We can convert these to a surface charge density by dividing each value 
by the billboard distance represented by the line charge. The first line 

1 2 24 m 

charge serves approximately 3/2 ( — j-g ); the last line charge 

serves 1/2 - ); and all others are associated with a length ~ )• 

2 

Figure 34 is a plot of charge/unit area (u Coul./m ) on the billboard as a 
function of length (northward) along the billboard surface. 

When an additional line charge in placed at the position of the lightning 
static, and all of line charge values are adjusted to the new configuration, 
we find the simulated potential function around a protected billboard - Figure 
35. The placement of the static in this example is based upon the discussion 
in Section 1.2.3. , with L =0.98m, corresponding to a = 40°. The charge/unit 
length for the static is 4.6 u Coul ./m. The charge/unit lengths for the ten 
billboard line charges in (u Coul./m) are: 

0.315, 0.47, 0.51, 0.57, 0.87, 0.89, 0.90, 1.35, 1.78, 2.1. 

These line charges may be compared with the unprotected billboard charges 
corresponding to the solutions of Figure 35. The protected billboard charges 
approach approximately one-half of the corresponding unprotected charges. 

The line charges used to simulate the rectenna are normalized to a 
charge/unit area through division by the associated lengths, as previously 
described, to obtain the induced charge distribution on the protected rectenna 
bill board. 

2 

Figure 36 is a plot of charge/unit area in uCoul./m as a function of the 
distance (northward) along the billboard face. 
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FIGURE 33 
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VIII. COMPUTATION OF LIGHTNING ELECTRIC FIELDS 


In section VII, a rectenna was simulated in the presence of a uniform 
electric field of 100,000 Volts. The induced surface charges derived from the 
simulation are directly proportioned to the imposed electric field strength. 

In this section we describe a computer program that was written to derive 
values for the lightning-produced electric fields as a function of time and of 
distance from "ground zero" - the point of strike. We have run the program 
for a range of lightning parameters obtained from actual measurements reported 
in the literature. 

The program computes the contribution to the electric field from the 
thundercloud charge center participating in the cloud-to-ground flash, the 
charge on the lightning channel, and the images of these charges. All charges 
are allowed to vary with time in a manner consistent with observations [Ter- 
restial Environment (Cl imatic) Criteria Guide! ines for Use in Aerospace 
Vehicle Development, 1977 Revision ; Edited by John W. Kaufman, NASA Technical 
Memorandum 78118]. 

Figure 37 displays the relevant equations and conf i gurati ons covering the 
leader phases of the computation. 

In Figure 38 the equations and conditions during the return stroke por- 
tion are shown. The program used in computing the fields is provided in the 
appendix. 

The material following Figure 38 provides the tabular and graphic data 
used in these computations for the return stroke phase. These data are 
contained in Figures (39-44) inclusive. 

The output of the computer program is a "blow-by-blow" history of the 
electrical field at a specified distance from ground zero as a function of 
time. Figure 45 displays one section of the output from one of the computer 
runs. This corresponds to a worst-case situation, 10 meters away from the 
very- severe-model . The units of time are seconds(along the abscissa), and the 
units of the ordinate are kilovolts per meter. 

Table 8.4 in figure 46 provides a summary of the output for the various 
computer runs. Listed are the peak negative fields, the peak positive fields 
(when positive fields occur), and the £ and AT for the portion of the flash 
with the peak rate of change of electric field. 

These values are our input data to the computation of diode failure when 
used in conjunction with the induced surface charge results of the rectenna 
electrostatic simulations. 
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INITIAL SPECIFICATIONS 
TFliPORAL FUNCTIONS: 

X = Y 0 - V l t 
q = Q 0 - P L (Y - X) 


'o 5 km) 

Q 0 (~-10 Coui) 

V L (^10^ m/s) 

t = 0 ) Y = Y 0 
Q l (~-5 Cou 1 ) 

„p-p l -«A 


SOLVE FOR E l (t,d) FOR t£ t l WHERE 

t l - (Y 0 -X L )/ V L 

X L (~50 METERS) 

FOR t>t l ) E l (t,d) = E^ (t l ,d) 



/ f / 


///////// 
■1 


E _ ) (D 2 + l 2 ) 1 ' 1 (D 2 + Y 2 ) 1/2 


1/2 


2 QY, 


'ITT£ 0 (D 2 + Y 0 2 ) 


3/2 


FIGURE 37 




' ”"l | for y ?y 0 ana p?-p l 

Q = pdT+ Q L J ^ q < _ ( q o _ Q^) 

SOLVE FOR E r (t',b) FOR t ? t l OR t' ? o 

TOTAL FIELD E t (t,d) = E L (t l ,d) + E R (t,d) 

TERMINATE COMPUTATION WHEN Q ? -(Q 0 - Q L ) 



///////// ///////// 


E = 2 £_ 1 

TttSd / 

•2'PYn ■ 

i nr£ 0 (D 2 + y 0 2 ) 3/2 



FIGURE 38 



DETAILS OF A VERY SEVERE LIGHTNING MODEL (MODEL l) 


Stage 

Key Points 

Rate of Current Change 

Charge Passing 

1. First Return 
Stroke Surge 

t = 0 i = 0 

t = 2 ps i = 200 kA 

t= 100 ps i= 7 kA 

| Linear Rise - 100 kA//is 
| Linear Fall - 193 kA in 98 ps 

0.2 C* 
~ 10. 2 C 

2. First Stroke 
Intermediate 
Current 

t = 100 ps i = 7 kA 

t = 5 ms i = 1 kA 

| Linear Fall - 6 kA in 4. 9 ms 

19.6 C 

3. Continuing 
Current — 
First Phase 

t = 5 ms i = 1 kA 

t = 55 ms i = 400 A 

| Linear Fall - 600 A in 50 ms 

35.0 C 

4. Continuing 
Current — 
Second Phase 

t = 55 ms i = 400 A 

t = 355 ms i = 400 A 

Steady Current 

120.0 C 

5. Second Return 
Stroke Surge 

t = 355 ms i = 400 A 
t = 355. 002 ms i = 100 kA 
t = 355. 1 ms i = 3. 5 kA 

} Linear Rise ~ 50 kA/ jus 
| Linear Fall - 96. 5 kA in 98 n s 

i 

~ 0.1 c 

~ 5.1 C 

6. Second Stroke 
Intermediate 
Current • 

t = 355. 1 ms i = 3. 5 kA 

t = 360 ms i = 500 A 

| Linear Fall - 3 kA in 4. 9 ms 

9. 8 C 


* Coulomb (C) is the quantity of electricity transported in one second by a current of one ampere. 
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CURRENT INCREASING 



DIAGRAMMATIC REPRESENTATION OF A VERY SEVERE LIGHTNING MODEL 
(MODEL l) (Note that the diagram is not to scale) 
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DETAILS OF A 98 PERCENTILE PEAK CURRENT LIGHTNING MODEL (MODEL 2) 


Stage 

Key Points 

Rate of Current Change 

Charge Passing 

1, First Return 
Stroke Surge 

t = 0 i = 0 

t = 5 pis i 100 kA 

t = 100 pis i 3. 5 kA 

1 

) 

) 

1 

Linear Rise - 20 kA/pis 
Linear Fall - 96. 5 kA in 95 pis 

0.3 C 
~ 4.9 C 

2. First Stroke 
Intermediate 
Current 

t = 100 pis i 3. 5 kA 

t = 5 ms i 500 A 

' 

, 

i Linear Fall - 3 kA in 4. 9 ms 

9. 8 C 

3. Continuing 
Current — 
First Phase 

t = 5 ms i 500 A 

t = 55 ms i 200 A 

1 

) 

Linear Fall - 300 A in 50 ms 

17.5 C 

4. Continuing 
Current — 
Second Phase 

t = 55 ms i 200 A 

t = 355 ms i 200 A 

Steady Current 

60 C 


FIGURE 41 




CURRENT INCREASING 



, DIAGRAMMATIC REPRESENTATION OF A 98 PERCENTILE PEAK CURRENT 
LIGHTNING MODEL (MODEL 2) (Note that the diagram is not to scale.) 
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DETAILS OF AN AVERAGE LIGHTNING MODEL (MODEL 3) 


Stage 

Key Points 

Rate of Current Change 

Charge Passing 

1. First Return 
Stroke Surge 

t = 0 i = 0 

t = 5 /us i = 20 kA 

t = 100 ms i = 2 kA 

| Linear Rise - 4 kA/jas 
| Linear Fall - 18 kA in 95 /is 

0.1 C 
~ 1.0 C 

2. First Stroke 
Intermediate 
Current 

t = 100 n s i = 2 kA 

t = 5 ms i = 300 A 

| Linear Fall - 1. 7 kA in 4. 9 ms 

5.6 C 

3. Continuing 
Current — 
First Phase 

t = 5 ms i = 300 A 

t = 55 ms i = 100 A 

| Linear Fall - 200 A in 50 ms 

10.0 C 

4. Continuing 
Current -- 
Second Phase 

t = 55 ms i = 100 A 

t = 355 ms i = 100 A 

Steady Current 

30.0 C 


FIGURE 43 
























CURRENT INCREASING 



DIAGRAMMATIC REPRESENTATION OF AN AVERAGE LIGHTNING MODEL 
(MODEL 3) (Note that the diagram is not to scale.) 


FIGURE 44 
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TABLE 8,4 


VERY SEVERE MODEL 


98 PERCENTILE MODEL 


AVERAGE MODEL 


Distance 

Peak 

Negative 

Peak 

Positive 

AE/AT 

Peak 

Peak 

Negative 

Peak 

Positive 

AE/AT 

Peak 

Peak 

Negative 

Peak 

Positive 

AE/AT 

Peak 

10 in 

-8. 5X10 5 

2.8X10 6 

2. 2X10 6 

-5.95X10 5 

1.81X10 6 

6. 46X10 5 

-5. 09X10 5 

1.30X10 6 

5. 68X10 5 



1.2X10' 5 



3.00X10" 6 



2. 59X10" 5 

50 m 

-5. 7X10 5 

1.7X10 5 

4.37X10 5 

-3. 88X10 5 

1.04X10 5 

3. 59X10 5 

-3. 10X10 5 

6. 1X10 4 

1.14X10 5 



2. 2X10' 5 



2. 5X10" 5 



2. 50X10" 5 

100 m 

-3. 49X10 5 

2.49X10 4 

2. 15X10 5 

-2. 36X10 5 

N/A 

1.75X10 5 

-1.85X10 5 

N/A 

5.47X10 4 



2. 2X10' 5 


2.5 X10' 5 


3.5X10" 5 

500 m 

-8.94X10 4 

N/A 

3.79X10 4 

-6.15X10 4 

N/A 

2.96X10 4 

-5.12X10 4 

N/A 

N/A 


3.2 X10* 5 


4. 5 X 10’ 5 


1000 M 

-5.35X10 4 

N/A 

1.69X10 4 
4. 2X10' 5 

-2.61X10 4 

N/A 

N/A 

-3.29X10 4 

N/A 

N/A 
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IX. COMPUTATIONS OF DIODE FAILURE 


We are now to the point of having generated all of the data that are 
required to evaluate the conditions under which the microwave rectifier diodes 
will fail due to induced currents from nearby lightning flashes. For a 
given £ and AT (from Table 8.4) we obtain from Figure 31 the power required 
for diode failure and from Figure 32 the induced charge/unit area on the 
rectenna surface. We assume that a diode designed to operate at 67 V will 
have a breakdown voltage of about 100 Volts. 


The surface area of the rectenna that has an induced surface charge of 
the size sufficient to cause diode failure is then computed from comparison 
with areas of the rectenna served by individual diodes and by series strings 
of diodes. Sample computations follow. 


SAMPLE COMPUTATION OF DIODE FAILURE 
(98TH PERCENTILE - 10 METER - NO PROTECTION) 

1. 98 percentile model - 10 meters: AT = 3 x 10" 6 and AE = 6.46 x 10 5 . 

2. Expected diode failure power from Figure 30: 250 Watts. 

3. Energy dissipated in the diode: 250 Watts x 3 x 1"° s = 7.5 x 10" 4 

Joul es. 

4. Charge transferred across 100 Volts diode breakdown voltage = 7.5 x 10" 6 

Coul ombs. 

5. From AE ip step 1 and figure 37, the induced charge/unit area 

= 3 xlO" c/nn x 6.46 = 19.38 x 10"° c/m. 


6. From steps 4 and 5, the rectenna area with surface charoe equivalent to 

the charge required to cause diode failure is: 0.39 m . 

7. Area served by diodes: rectenna center, . 


25 watts 


= 0.11 m ; rectenna edge. 


25 watts 


2.5 m 


o v* x± m 5 i cLLCima cuyc , p 

230 w/m 10 w/m 

8. Compare 6 with 7: single diode configuration near rectenna center is 

safe. Single diode conf iguration near rectenna edge is vulnerable. 

9. However , the diodes are to be put in series (597 to a string) hence the 

diodes near the bottom must carry all of the induced current to the 
entire string. For these bottom-string diodes the area served with 
respect to the induced charge is: rectenna center, 60 m , rectenna edge, 

1400 m 2 . 

10. To protect against the 98 percentile flash within 10 meters of ground 

zero would require fast surge protection diodes (back to back zeners) on 
all diodes in the rectenna. This extent of protection may not be cost 
effective; however the considerations in Section X indicate that simpler 
protection arrangements will probably be effective near the rectenna 
center. 


FAILURES PRODUCED BY THE AVERAGE LIGHTNING FLASH 


The situation considered here is the extent of the protection required 
for an "average" lightning flash if we are willing to accept losses from the 
extreme cases. 
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The computation sequence follows the same procedure described immediately 
above. Here we use data for the average flash from Table 8.4 at a 10 m 
distance from ground zero. 


SAMPLE COMPUTATION OF DIODE FAILURE 
(AVERAGE FLASH, 10 M, WITH “STATIC" PROTECTION) 

1. From Table 8.4: ZE = 5.68 x 10^ v/m; AT = 2.59 x 10" ^s. 

2. From Figure 6.1: _80 watts. _ ^ 

3. 80 w x 2.59 x 10 s - 2 x 10" Joules. 

4. 2 x 10 -b coulombs. _ r ? 

5. From 1 and Figure 38: 1.5 x 10" D x 5.68 = 8.52 x 10" coul/m . 

6. From 4 and 5: Area = 2.35 m . 

7. Since the rec£enna area served by individual diodes even on the 

edge < 2.5 m , the individual diodes are self-protecting and able to take 
an "average" lightning flash. 

8. However , when arranged in a series stack of 597, the diodes at the 
bottom of the stack must conduct the induced currents for the whole 

stack. The diodes cannot safely carry these currents. 


60 



X. LIGHTNING PROTECTION FOR SERIES DIODE STRINGS 


As demonstrated in Section IX, the connection of microwave rectifier 
diodes in series requires special lightning protection considerations. We 
cannot make specific recommendations for these protection devices at this time 
because the rectenna current design is not advanced to the point that allows 
such detailed analysis. Rockwell International has provided us with an equiv- 
alent circuit for the rectenna; a slightly modified form of that circuit is 
shown in Figure 46. We have assumed that the series connections are to be 
made at the points indicated by the large spots and that the output filler 
operates around 30 Hz. A series string of rectenna elements of this design 
can be protected with a variety of methods. One cost-effective means is a 
spark gap arrangement incorporated in the diode feedthroughs, or the output 
filter inductors, or on the billboard configuration itself. 



RECTENNA EQUIVALENT CIRCUIT AT 2,45 GHz 



FIGURE 46 


XI. CLOUD-TO-GROUND LIGHTNING DISTRIBUTION IN THE UNITED STATES 


In order to have a working estimate of the hazard presented by lightning 
to rectennas, we need to know the cl oud-to- ground lightning flash density for 
various possible rectenna sites in the United States. The cloud-to-ground 
lightning flash density (in #/knr for example) is not a parameter that is. 
measured as a climatological variable. We have found it necessary to use the 
number-of-thunderstorm days as a proxy variable because it is available as a 
climatological variable. Figure 47 gives contours of annual number-of- 
thunderstorm days. 

XI. 1. Pierce Conversion Formula 

Several attempts have been made to derive a conversion formula to convert 
thunderstorm days into the flash density by using lightning flash counters in 
research areas for correlation with the count of thunderstorm days. The best 
of the various conversion formulas is that due to E.T. Pierce ("A Relationship 
Between Thunderstorm Days and Lightning Flash Density," Trans. AGU, 49 , 686, 
1967.) The Pierce formula (as does most others) has a quadratic term, which 
reflects the relationship between frequencies of local storms and storm 
intensity. In addition, the formula utilizes the monthly thunderstorm days as 
opposed to the annual average in order to incorporate seasonal effects in the 
conversion formula. 


This formula is 

2 2 4 

= aT M + a T M » 

where: T M = monthly number of thunderstorm days and is the monthly ground 

flash density (#km^/Mt.) The parameter a is, 

a = 3 x 10" 2 

If a is the annual ground flash density (# knTVyr.), then 

W=1 M 

XI. 2. Climatological Data -- Number of Thunderstorm Days 


The inputs needed to compute the U.S. Distribution of ground lightning 
flash density are: (1) The monthly number of thunderstorm days for all U.S. 

stations recording these observations, (2) the coordinates of the observing 
sites, and (3) the computer software to compute the density and display the 
results geographically. 

Items 1 and 2 were obtained from "Local Climatological Data - Annual 
Summaries for 1977" published by The National Oceanic and Atmospheric Adminis- 
tration on magnetic tape. The geographic plotting software of Item 3 was 
obtained from The National Technical Information Service, and the computer 
programming was done by J.L. Bohannon at Rice. 
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A detailed list of flash density for all of the stations used is provided 
in the Appendix. 

Note the hot spots on the contours in Figure 48 that result when stations 
are located near geographic features that promote local thunderstorms. There 
are probably other similar hot spots in the U.S. that do not show up on this 
display because of the absence of an observing station nearby. 
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APPENDICES 


Computer programs developed under this contract. 

All programs are in FORTRAN H, unless otherwise specified. All of the 
programs were run on an IBM 370/155 and/or an Itel AS/6 computer. 
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Apendix A_ 


Computer Program PANEL: 

A Computer Model of the SPS Plasma Interaction 


The following pages are the listing of the program "PANEL," written to 
model the interaction of a high voltage solar array with an ambient Maxwellian 
plasma. The program was originally written by Dr. Lee W. Parker and was 
modified for application to the SPS problem by David L. Cooke. 
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79 

PANEL 

N 1 r 

4 

END 

JOB 

24 04 

B IN 

7 

COOKE 

437075SCPS 

4:02:00 

PM 

27 

MAR 

79 

PANEL 

N t E N 


END 

JOB 

2404 

BIN 

7 

COOKE 

437075SCPS 

4:02:00 

PM 

27 

MAR 

79 

PANEL 

x J ER3 

A 

END 

JOB 

2404 

BIN 

7 

COOKE 

437075SCPS 

4:02:00 

PM 

27 

MAR 

79 

PANEL 
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> 

w 


ISN 

1SN 


TSN 

ISN 

ISN 

ISN 
ISN 
T SN 
I SN 
ISN 
TSN 
ISN 
I SN 
ISN 
ISN 


I SN 
ISN 
I SN 
TSN 
I SN 
1 SN 
TSN 
ISN 
1 SN 

I SN 
1 SN 
ISN 
1 SN 
I SN 
TSN 
ISN 
ISN 
1 SN 
TSN 
I SN 
I S N 
I SN 
ISN 


I SN 
1 SN 
1 Sn 
1 SN 
l Sn 
I SN 
1 SN 


0002 

0003 


0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 


0017 

0010 

001 9 
0020 

002 1 
0022 

0023 

0024 

002 5 

0026 

0027 

0020 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

003 7 
00 3 0 
0039 


0040 
0042 
004 3 

0044 

0045 

0046 

0047 


COMPILER OPTIONS 
C 


NAME = HAIN 1 OPT=O2,LINECNT = 60.SIZE=OOO0K, 

SOURCE, EBCDIC, NOLIST, NOD ECIC.LD AD, MAP ,NOEOIT,NOID,NOXREF 


C SOLAR PANEL PROBLEM 

CQMMQN/CP/NPRINT,NPrS,MA,MB»ME»KMAX»XPT»YPT*ALlfBEl»EV»SMACH» 
1 TVOLTStCUR, XMETER 

COMMON/ OK/IIM, IIP, JJM,JJP*KK«NTOT ,IV»JV»II *JJ,M*N,VP(30), 
1XYZC2080,3),VVC30,20,10),XPC30),XMC10),YP{20},YM(10),Z2C10), 
2XX(4Q3»YY(309* ILX , IU X , KUK , HOC , MB 13 , V RF , NF P S , SKPRF L , SKPLST 
C0MMDN/FLD/XC2080.2) , COE FC 2000, 7 ) , IN0XC2 0 80, 6 ), SKPCO 
C0MM0N/C0/PV0LTS,XMACH|0ENSI,NN,PARTCL(2),PARTlC29,PART2C2) 
COMMON/INTER/ IN T ,1 1 A, JJA, KKA, IGQU.T, JG 0 UT , K GO UT, X A , Y A , Z A , 

IX I C 309 ,YJC2Q)«ZKC 10) 

OIMENS ION D A T E C 2 0 ) 

DIMENSION VFCC4),IFC4),JFC4),KFC4) 

INTEGER SXPRFL, SKPLST, SXPCO 
NF(IX,JX,KX)=IXfII*CJX“l)+II*JJ*CKX“l) 


L- 5 
M= 6 


100 READCL, 9999, END=99) DATE 
9999 FORMAT C 2 0 A 4 ) 

2 URITECM, 9998) DATE 

9998 FORMAT (42H1S0LAR PANEL ELECTRIC FIELD AND CURRENTS. ,20A4) 

C READ GEOMETRIC PARAMETERS 

C 

RE AD CL ,111) IIP.IIM, JJP, JJM,KK,IV, JV 

I I — 1 1 M + T IP-1 
JJ= JJM+JJP-1 
NTOT=II*JJ*KK 

READCL. 222) C X P C I ) , I = 1 , 1 1 P ) 

READCL ,222) C X MC 13 , 1=1 , I I M) 

READCL.222) ( Y P C J ) , J = 1 , J J P ) 

READCL, 222) CYMcJ) , J=1 , J JM) 

R E A 0 C L , 2 2 2 ) CZZCK),K=1,KK) 

C READ PANEL POTENTIALS 

RI':A0CL,116)CVPCI),I = 1.IV),VRF 

READCL , I U)SKPRFL,SKPLST , IL X , IUX , KLK , K UK , MBC ,MBD , NFP S , SKP CO 
DO 140 NPt=l,NTOT 
XCNPC, 1)=0 
140 XCNPC, 2)=0 

J I M 1 — I1MI-IV-1 
JJM1 = J JM ► J V“1 
DO 150 I = IIM.IIH1 
DO 150 J = JJM.JJMl 

III = jUl-IIM 
N = NFCI ,J,1) 

X C N , 1 ) = VPCIII) 

XCN.2) = 1 
150 CONTINUE 
C 

C CONSTRUCT REFLECTORS 

IFCSKPRFL.EO.DGO TO 163 
DO 160 I = ILX, IUX 
DO 160 K = KLK.KUK 
JW = MBC-K 
NW = NFC I, JW, K) 

X CNW, 1 ) = VRF 
X ( NW, 2 ) = 1 



> 


ISN 

0048 


ISN 

0049 


ISN 

0050 


I SN 

005 1 


ISN 

0052 

160 

ISN 

0053 


ISN 

0054 

231 

C 

C 

ISN 

0055 

163 

ISN 

00 57 


ISN 

005 8 

118 

ISN 

0059 


ISN 

0060 


ISN 

0061 

119 

T SN 

0062 


ISN 

0063 

117 

ISN 

0064 

165 

ISN 

0065 


ISN 

0066 


ISN 

0067 


ISN 

0068 

170 

ISN 

0069 

220 

I SN 

0070 


T SN 

007 1 


ISN 

0072 


ISN 

0073 


I SN 

0074 


ISN 

0075 


T SN 

0076 


ISN 

0077 


ISN 

0078 


ISN 

0079 


ISN 

0080 


ISN 

000 1 

C 

I SN 

0082 

111 

ISN 

. 0083 

113 


ISN 

0084 

116 

ISN 

0085 

222 

ISN 

0086 

223 

ISN 

0087 

224 

ISN 

008 8 

225 

I SN 

0089 

226 

ISN 

0090 

227 

T SN 

009 1 

228 

1 SN 

0092 

229 

T SN 

0093 

230 

ISN 

0094 

240 

I SN 

0095 

241 


JW = K+M8D 
NW = NF(I,JW,K) 

X ( N W , 2 ) = 1 
XCNW.l) = VRF 
CONTINUE 

WRITECM, 231) VRF 

FORMATC//1X, 'REFLECTOR POTENTIAL = '.1PE15.5) 

READ ADDITIONAL FIXED POTENTIALS 
IF(NFPS.LE.3)GO TO 220 
WRI TEC M, 118) 

FORHATC//' ADDITIONAL FIXED POTENTIALS'/ 

14 (6X,* POT* ,7X,' I* ,3X . *J' ,3X, *K» )) 

DO 170 NOC = 1,NFPS,4 

READ(L»119)(VFC(I),IF(I)»JFCD»KF(I), 1=1,4) 

FOR HA T (4 (EB. 0*314)) 

WRITECM, li7)<VFC(I),IFCI),JF(I),KF(I>, 1 = 1,4 ) 
F0RMATC/4C3X, IPEIO.2,314)) 

DO 170 1=1.4 

NN = NF(IFCI).JFCI).KFCI)) 

XCNN, 1 ) = VFC ( I ) 

X(NN«2)=1 
CONTINUE 
CONTINUE 
IVP=I V+l 
JVP=JV+1 


WRITECM, 113)IIP,IIH, JJP, JJH.KK, IV, JV 
WRITECM, 223) CI,XPCI),I = 1, IV) 

WRITECM, 224) CI,XPCr),I=lVP,IIP) 

WRITECM.22S) CI.XM(I),I = 1* IIM) 

WRITECM, 226) CJ, YP(J), J=l, JV) 

WRITE CM, 22 7) C J, YPC J) , J = JVP , JJP) 

WRITECM, 228) ( J , YM ( J) , J=1 , J JM) 

WRITECM, 229) ( K , ZZ ( K) , K.= 1 , KK) 

WRITECM, 230) CXPCI), 1=1 , IV) 

WRITE CM, 24 1)C VP Cl), 1 = 1, IV) 

FORMAT Cl 61 5) 

FOR MAT C//1X, I 3, 1 8H POSITIVE X-VALUES/ 

1 1X.I3.18H NEGATIVE X-VALUES/ 

2 IX, 13, 18H POSITIVE Y-VALUES/ 

3 IX, 13, 18 H NEGATIVE Y-VALUES/ 

4 IX, 13, 25 H Z-VALUES C P 0 S I T I VE ONL Y ) / 

5 IX, 13,33 H POSITIVE X-VALUES DEFINING PANEL/ 

6 IX, I3.33H POSITIVE Y-VALUES DEFINING PANEL) 

FORMAT C8 El 0. 0 ) 

FORMAI ( 1 6E5. 0) 

FORMAT C//1X.27HX-VALUES POSITIVE ON P A NE L = / C 1 3, 1 PE 15 . 4 )) 

FORMAT C//1X, 35HX-VALUES POSITIVE OUTSIDE OF PANEL=/CI 3, 1PE15.4)) 
FORMATC // IX, 18HX- VALUES NEGATIVE=/CI3,1PE15.4)) 

FORMAT C//1X.27HY-VALUES POSITIVE ON P AN E L = / C 13, 1 PE 1 5 . 4 ) ) 

FORMAT C/ /1X,35HY-VALUES POSITIVE OUTSIDE OF PANE L = / C 1 3 , 1 P E 15 . 4) ) 
FORMAT C//1X, 18HY-VALUES NEGATIVE = /CI3,1PE15.4)) 

FORMAT C//1X, 37 HZ -VALUES (POSITIVE ONLY) ABOVE P ANEL = /C 13, 1PE15.4)) 
FORMAT C////1X .25HARRAY OF PANEL POTENTIALS// 

1 15X.3HX =,3X,(8(F8.4»4X)/20X)) 

FORMAT C/1X.2HYC, 12, 2H)=,F8.4,6X,(8( 1PE 12.4 )/20X)) 

FORMAT C 3X, 'ALL Y,',5X,(8(1PE12.4)/20X)) 



PAGE 003 


1 • 1515) 

/ix,36HNPRINT,NPTS,MA,MB,ME,KMAX .PROBNO =,616.110/ 
NUMBER =, F9.1, 9X, 13 H TEMPERATURE =, F9.1, 6H VOLTS, 9X 
F9.1, 7H PER CC, 9X, 6HMASS =, F9.0.11H ELECTRONS 
II SCALE =,F9. 1, 30H METERS = X-DIMENSION OF PANEL) 

? JUS INGLE SPACE POINT. X =,F10.S,5X, 3HY =,F10.5) 

3 5 H SINGLE ENERGY ( MONOENE RGE T IC ) . E =,F10.5, 6H VOLTS) 

2 T IS INGLE TRAJECTORY. X =,F10.5,5X, 3HY =,F10.5/ 

1.1 ALPHA = » F2 0 . 8 , 8H DEGREES/ 

LI. BETA =, F20. 8 , 8H OEGREES/ 

R GY =»F20.5 , 6H VOLTS) 

T2HRANOOM THERMAL CURRENT DENSITY =,1PE13.6, 

FI SUUARE METER, F0R.2A5) 

II INTERFACE X-VALUES/(I3,1PE15.6)) 

I UHINTERFACE Y-VALUES/(I3,1PE15 .4 )) 

LUlt INTERFACE Z-VALUES/CI3, 1PE15.4)) 

!' "22 H -- CURRENTS AND POWER)) 


, 

/ 


[ I 


i,0 TO 620 
I 


> 

cn 


0 TO 680 


VC/)- YnfO - "5 

i io 0 


/ViT 


o 


o 


1 “1 ) *XX(I )) 


J~1 )♦ YY( J)) 


ISN 

0141 


ZKC1)=ZZC1) 

ISN 

0142 


ZICCKKA) = ZZCKK) 

T SN 

0143 


DO 560 K=2,KK 
ZKCK)=.5*CZZCK-l) + ZZCIO) 

1 SN 

0144 

560 

1 SN 

014 5 


WRIlECMj 561) CI,XICI),I=1,IIA) 

ISN 

0146 


URITECM.562) < J» Y JC J) , J=1 , JJA) 

I SN 

0147 

r 

URITECM.563) CK , ZK OO , K= 1 , KK A 3 

ISN 

0148 

L» 

DO 600 N=1,NT0T 

ISN 

0149 


CALL FINDCIFIND, JFINO.KFIND) 

ISN 

0150 


XYZCN, 1)=XXCIFIN0) 

ISN 

0151 


XYZCN, 2)=YYCJFIND) 

ISN 

0152 


XYZCN, 3)=ZZC KFINO) 

ISN 

0153 

600 

r 

CONTINUE 

ISN 

0154 


IFCSKPLST.EQ.l) GO TO 660 

ISN 

0156 


NFP P = C NTOT/3003+1 

ISN 

0157 


DO 650 I P=1 • NFPP 
UR I IE ( M, 900 5 ) 

ISN 

0158 


ISN 

0159 

9000 

F0RMATC1H1/6X,1HN.3X,4HXCN),2X,4HY(N),2X,4HZ(N)//) 

ISN 

0160 


CALL LISTC2.IP) 

ISN 

0161 

650 

CONTINUE 

ISN 

0162 

660 

r 

CONTINUE 

ISN 

0163 

l* 

DO 700 J=1,JJ 

T SN 

0164 


DO 700 1=1,11 

ISN 

0165 


N = NFCI.J.l) 

I SN 

0166 


vvci.j.d = xcN.n 

ISN 

0167 

700 

r 

CONTINUE 

ISN 

0160 

L 

K= 1 

> ISN 

0169 


URITE(H,8000) K.ZZOO.CXXCI), 1=1,119 

5a ISN 

0170 


DO 750 J=1,JJ 

ISN 

0171 


UR I TECH, 24 0) J, Y Y ( J) , ( VVC I , J , K) , 1=1 , I I) 

ISN 

0172 

750 

r 

CONTINUE 

ISN 

0173 

L 

r 

CALL FIELO 

ISN 

0174 

u 

00 800 K=1,KK 

T SN 

0175 


DO 800 J=1,JJ 

ISN 

0176 


DO 800 1=1,11 

ISN 

0177 


N=NFCI , J , K) 

ISN 

0178 


VVC I , J » K.) = X ( N , 1 ) 

ISN 

0179 

800 

r 

CONTINUE 

ISN 

0180 

L 

DO 900 K=l,KK 

ISN 

0181 


URITECH.8000) K,ZZ(IC),CXX(I),I = 1 .11) 

ISN 

0182 

8000 

FORMAK 26H1ARRA Y OF POTENTIALS AT Z C , I 2 , 2H)= , F8 . 4// 

L 15X.3HX =,3X,C8CF8.4,4X)/20X)) 

ISN 

0183 


DO 850 J=1,JJ 

ISN 

0184 


UR I TE C H, 240) J»YYCJ),CVVCI»J»IO»I = l»II) 

I SN 

0 1 B 5 

850 

CONTINUE 

ISN 

018 6 

900 

r 

CONTINUE 

ISN 

018 7 

L 

NP R 0B= 0 

ISN 

0188 

1000 

RE AOCL ,333,END = 99) N PR I NT , NPTS , M A , MB , ME . K HAX .MORE 
RE ADC L , 1 16) SMACH, T VOL TS .DENCC.X MASS, XHE TER 

ISN 

0189 

1001 

ISN 

0190 


NPROB = NPROB<- 1 



I SN 0191 
ISN 0192 

1 SN 0193 
ISN 019!) 
ISN 0197 
1 S N 019 9 
ISN 0201 
ISN 0203 
ISN 0205 
ISN 0207 
ISN 0209 
ISN 0211 
ISN 0213 
ISN 0214 
ISN 0215 
ISN 0217 
ISN 0218 
ISN 0219 


WRITECM.999) 

URITEC Mi 444 0 N PRI N T, NPT S , MA , MB , HE , KM AX , N PR OB , SMACH .TVDLTS .DENCC. 
1 XNASS iXMETER 

IFCNPTS.EQ.O.OR.ME.EQ.O.OR.MA.EQ.O) READCL»222)XPT,YPT,AL1,BE1,EV 
IFCNPrS.EQ.03 UR 1 TEC M, 445 1 XPT.YPT 
IFCME.EQ.03 WRITECM,446) EV 

IFCMA.EQ.O) URITECM,447) XP T , Y PT , AL 1 , BE 1 , E V 
IFCMA.GT .O.AND.XMASS.LE.O.) STOP 

IF(MA.GT.O) CUR=2.68£-8*DENCC*SQRTCA8SCTVQLTS)/XHASS) 
IFCTVOLTS.GT.O.3 P AR TCL C 1 ) = P ART 1 ( 1 ) 

IFCTVOLTS.GT .0.3 PARTCLC2)=PART1C2) 

IFCTVOLTS.LT.O.) PARTCLC1)=PART2C1) 

IFCTVOLTS.LT .0.1 PARTCLC2 )=PART2C2) 

URI TEC Mi 448) CUR.PARTCL 
CALL POWER 

IF CMORE.GT.O) GO TO 1000 
GO TO 100 
99 STOP 
END 
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DATE 


I SN 

0002 

ISN 

0003 

ISN 

0004 

ISN 

0005 

ISN 

0006 

ISN 

0007 

ISN 

0008 

I SN 

0009 

1 SN 

0011 

I SN 

0013 

ISN 

0014 

ISN 

0015 

ISN 

0016 

T SN 

0017 

ISN 

0018 

ISN 

0019 

ISN 

0020 

ISN 

0021 

ISN 

0022 

I SN 

0023 

ISN 

0024 

ISN 

0025 

ISN 

0026 

ISN 

0027 

I SN 

002 9 

ISN 

0030 

ISN 

003 1 

ISN 

0032 

ISN 

0033 

ISN 

0034 

ISN 

0036 

ISN 

0037 

ISN 

0038 

I SN 

0039 

I SN 

0040 

I SN 

0042 

ISN 

0044 

ISN 

0046 


COMPILER 


r, 

c 

c 


999 

C 


OPTIONS - N A ME = M AI N , 0 PT = 02 , L IN E C N T = 60 , S 1 1 E = 0 OOOK , 

SOURCE, EBCDIC, NOLIST.NO DECK, LOAO, M AP .NOEDIT, NO 10, NOXREF 
SUBROUTINE ORBIT 

STEP ACROSS 3-0 BOX ASSUMING CONSTANT POTENTIAL WITHIN BOX 

COMMON/BK/IIH.IIP, JJH, JJP.KK.NTOT, IV, JV,II, J J, K, N . VP C30D , 
1XYZ(2080,3),VV(30,20,10D,XPC30),XM(10),YPC20D,YM(10D,ZZ(10D, 

2XXC 40 3 ,Y YC30 9 , ILX , IUX ,KUK , MUC.MBD ,VRF,NFPS , SKPRFL , SKPLST 
COMMON/ORB/XOOT, YDOT , ZDOT , X 1 , X2, Yl, Y2.Z1 , Z 2 , X , Y , Z . PH I , NT I HE . S A VE 
DIMENSION TIMEC6) « U( 3D ,UDOT C3),BC2, 3) 


T00M=3. 3333E +33 
ROUND = l.E-11 
TR-aWfB- - - 1 - v tr-Yr- 


(s Mo 


' s ) 


IF(XOOT.EO.O..AND. YDOT . E 0 . 0 . . AN 0 . ZDOT . EO . 0 . I WRITECM.999) 
IFCXDOT. EO.O..ANO.YOOT.EO.O.. AND. ZDOT. EO.O.) RETURN 
FORMAT ClX,38HSPEED-0 - HENCE PARTICLE DOES NOT MOVED 

ucn=x 

U C 2 D = Y 
U( 3D=Z 


UDOTC 1 D=XOQT 
U00TC2 D=YDOI 
UDOTC 3 3= ZDOT 

B(i,n=x i 

BC 2, 1)=X2 
B ( 1 , 2 D = Y 1 
BC2,2D=Y2 
B(1 ,3D=Z1 
B C 2 * 3D = Z2 




X 

5 


j 


*1 


100 


DO 101 N 2 = 1 , 3 

IF(U0QT(N2).EQ.0.D GO TO 101 
DO 100 N 1 = 1 , 2 
N R = N 1 + _2*CN2-1) 

TIME (NRD =TOGM 

TT= ( fl( N1 , N2 ) UCN2D D/U00TCN2D 
SS = UGN2D ' *- 31 0 0T(N2D*TT 

IF( SS CTiEvBt 1»N2D.AND.SS.LE.B(2,N2DD TIMECNRD=TT 
CONTINUE • 


101 CONTINUE 



E% 


ox- /( 

^ /i/J yU 


& 


O(H^) -t 


C FIND SHORTEST SIGNIFICANT TIME 
C 

T IM IN = TOQM 
DO 200 NR = 1 , 6 

IFCTIMECNRD.EO.TOOMD GO TO 200 / 

IFCTI MECNRD.GT.RO UNO. A ND . T I ME CNR D.LT . TIMIN D ,NTI ME =NR 

2 00 CONTINUE NR) * GT ‘ R0UNC> * ANlD * Ti:ME CNR;> * LT * TIMIn;) ( T IMI N = T I ME CNR D 
C ’ 

C ADVANCE TO APPROPRIATE END-POINT 
C 





. s 


. -pr 


ISN 

00 A 7 


X = X + XDOT*TIMIN X ^ Y 

x 04 r* 

T SN 

00 A 8 


’Y® Y + YOGI*! I M I N 


I SN 

004 9 

r 

/.= Z t ZDOr*T IMIN 


l SN 

0050 

Li 

XS AV=X 


ISN 

005 1 


YS A V = Y 


ISN 

0052 


ZSA V=Z 


ISN 

0053 

C 

IF C NTIME .EQ. 1 0 f=X 1\ r 


ISN 

0055 


IFCNTIME.EO.20 X-^XJU 


ISN 

0057 


IFCNTIME.EQ.30 Y = Y 1 


T SN 

0059 


IFCNTIME.EQ.40 Y=Y2 


ISN 

0061 


IFCNTIME.EQ. 50 Z=Z1 


ISN 

0063 

r 

ie.CUTIHE.EQ. 60 Z = Z 2 


ISN 

0065 


/' / DX=X-XSAV . 


T SN 

0066 

i 

f DY=Y-YSAV ) 

k^JOZ = Z-ZSAV^/ . 


ISN 

0067 

r 


I SN 

0068 

L 

IFCCNTIME.EO. 1.DR.NTIME.EQ.20. AND.ABSCDXO. GT.TROUNDO, 

NTIME=-1 

ISN 

0070 


IFCCNTTME.EQ.3.0R.NriME.EQ. AO. AND. A BSC DYO.GT.TR DUNDO/ 

NTIME=-2 

ISN 

0072 


IFCCNTIME.EQ.5.QR.NT I ME. EQ. 60. AND. ABSCDZO. GT.TROUNO) 

NTTME=-3 

ISN 

00 7 A 


TFCNTIME.EQ.-lO SAVE=XSA V \ 


ISN 

ISN 

0076 

0078 

r 

IFCN1IME .EQ. -20 SAVE=YSAV ^ 

IFCNTIME.EQ.-30 SAVE=ZSA V 


ISN 

0080 


RETURN 


ISN 

0001 


END 





> 
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OS/ 360 FORTRAN H 


ISN 0002 


ISN 0003 


ISN 0004 

ISN 0005 
ISN 0006 
ISN 0007 


COMPILER OPTIONS - NA ME-= M A I N, OP T= 02 . L I NE CNT=60 , SI ZE = 0000K, ., nvDC( - 

SOURCE, EBCDIC, N0LIST,N0DECK, LOAD , M AP ,N0EDIT , NO ID , NOXREF 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


0008 

0009 

0010 
0011 
0012 
0013 

0015 

0016 
0018 
0020 
0021 
0022 

0023 

0024 


SUBROUTINE OEN 

C ROUTINE FOR EVALUATING CURRENT-DENSIT Y INTEGRALS OVER VELOCITY SPACE 

C COMMQN/BK/IIM,IIP,JJM,JJP,KK,NT0T,IVfJV,IT,JJ»M»N»VP(30)» 

lXYZ(2080»3)»VV(30,20§10)»XP(30),XM(lQ),VP('p),YM(lp)»ZZC10), 
2XXC40 ), YYC30), ILX, IUX , KUK , MBC , MB D, VRF.NH‘5 , SKPRFL.SKPLST 
COMMON/CP/NPRINT,NPTS,MA,HB,ME,KMAX,XPT,Y.PI , AL1 ,BE1 * EV ,SMACH, 

I TVOLTS,CUft»XMEIER „ „ 

COMMON/CD/P VOLTS, X MACH , DENS T , NN, PARTCL C2) .PART 1C2),PART2C 2) 
C0MM0N/0RB/X00T,YD0T,ZDDT,X1,X2,Y1,Y2.Z1 , Z 2 • X , Y , Z, PH I, NT I ME . S AVE 
COMMON/INTER/ INF ,IIA, JJA.KKA, I GOUT, JGOUT .KGOUT.XA, YA.ZA, 1 

1X1(30), YJ(20),ZKC10) 

DIMENSION A(2),END1(2) ,END2(2),FATE(2) 

OATA END1/4HABS0,4HRBE0/,EN02/4HESCA,4HPES / 

XS AVE=XPT 
YS A VE = YPT 
TEMP= ABS(TVOLTS) 

IFCTEMP.LE.O.) WRITECM, 999) TEMP 

999 FORMAT C////1X, 38HTRQUBLE - NEGATIVE OR ZERO TEMPERATURE) 
IFCTEMP.LE.O.) RETURN 
IFCMA. EO.O.OR.ME.EO. 0) EE=EV/TEMP 
PI=3. 1415926536 
A(1)="1./SQRT(3.) 

A(2)="A( 1) 

M0STPS=0 
MSTEP=1000 


C 

C 


SET UP SUMS OVER TRAJECTORIES 
GO TO 250 


ISN 

0025 

C 

IFCMA. EQ. 

0) 1 

ISN 

0027 


JAMAX=2 


ISN 

0028 


JBMAX=2 


TSN 

0029 


KAN AX=MA 


ISN 

0030 


KB M AX = HB 


ISN 

0031 


NUMBER=MA 

*MB 

ISN 

0032 


IFCNN.EO. 

1) 

ISN 

0034 

990 

FORMAT C/1X, I 



r 

1 5HHENCE , 

14, 

ISN 

0035 

L 

IFCME.EO. 

0) 

ISN 

0037 


ME2 = 2* ME 


TSN 

0038 


JEM AX= 2 


ISN 

0039 


KE M AX = ME 


ISN 

0040 


IFCNN.EO. 

1) 

ISN 

0042 

988 

r 

FORMATC1X 

, 14 

ISN 

0043 

L 

GO TO 300 




C S 

INGLE VALUE 

OF 

ISN" 

0044 

200 

JEMAX=1 


ISN 

004 5 


K EM AX= 1 


ISN 

0046 


IFCNN.EO. 

1) 

ISN 

0048 

986 

FOR MAT C IX 

.31 


*4 

WRITECM, 990) 


GO TO 200 


MA, MB, NUMBER 


WRITECM 988) ME ME2 

,27H ENERGY INTERVALS AND HENCE, 14, 14H ENERGY VALUES) 


ENERGY 


WRITECM, 986) EV.EE 
H MONOENERGETIC CASE 
DIMENSIONLESS V A LUE , E 16. 4} 


WTIH ENERGY, 1PE16. 4, 30H VOLTS, OR 






.?• 


ISN 0049 


ISN 0050 
ISN 0051 
ISN 0052 
ISN 0053 
ISN 0054 
TSN 0055 
ISN 0056 
ISN 0057 
ISN 0058 
ISN 0059 


ISN 0060 
ISN 0061 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 


ISN 0072 
ISN 0073 
ISN 0074 
ISN 0075 

ISN 0077 
ISN 0078 
ISN 0079 

ISN 0080 
ISN 0081 


ISN 0082 
ISN 0083 
ISN 0084 
ISN 0085 
ISN 0086 

ISN 0087 

ISN 0089 


GO TO 300 
C 

C SINGLE TRAJECTORY ONLY 
C 

250 JAM AX= 1 
JBMAX=1 
JEMAX=1 
KAMAX=1 
KBMAX=1 
K E H AX = 1 
AL=AL1*PI/180. 

BE=BE1*P 1/180. ■ 

UR I TEC M, 984 ) AL1 , AL, BE 1 , BE , EV , E E 
984 F0RMAT(/1X,17HSINGLEI TRAJECTORY 

1/1X, 7HALPHA = « F20 . 8 ,12H DEGREES, 0R,F20.B , 8H RADIANS 
2/1X , 7HBEI A = , F2 0 . 8 ,12H DEGREES, 0R,F20.8 , 8H RADIANS 
3/ IX , 8HENERGY = , IP E16. 4 , 30H ' VOLT S, OR DIMENSIONLESS VALUE.E16.4) 
S IN A=S IN CAL) 

COSA=COSCAL) 

C 

C SUM OVER ENERGY, BETA, ANO ALPHA 
C 

300 CONTINUE 
DENST=0. 

DO 1001 KE = l .KEHAX 
DO 1001 JE = 1 , JE M AX 
DENS=0 . 

N0ESC=0 

DO 1000 KB=1 »K8MAX 
DO 1000 JB=1,JBMAX 
00 1000 KA=1 »KAMAX 
DO 1000 JA=1,JAMAX 

C INITIAL POSITION' 

C 

2 = 0 . 

X=XSAVE 
Y=YS AVE 

IFCMA.EQ.O) GO TO 320 
C 

CA=CACJA) + FL0ATC2*KA - 1 - MA33/FL0ATCMA3 
SINA=S0RTC.5*C1.+CA)) 

C0SA=S0RTC1. - SINA**2) 

C 

CBETA=CACJB) «• FLOAT C2*KB - 1 - HB33/FLQAT CMB} 

BE = PI*< 1 • *■ CBETA3 
C 
C 

320 XD0T=SINA*CQSC8E) 

YD0T=SINA*SINC8E) 

200T=C0S A 
INT =0 

CALL INTERP 
C 

IFCIGOUT.GE.l. ANO.IGOUT.LE.IIA. ANO. JGOUT . GE . 1 .AND. JGOUT.LE . J J A . 

1 ANO.KGOUT.GE.l.ANO.KGQUT.LE.KKA) GO TO 340 
330 WRITE CM, 99993 



I SN 

0090 

9999 

FORMAT (/////IX «*3H0NE OF THE IG-JG-KG INDICES IS OUT OF RANGE) 

ISN 

0091 


WRITE (M.B88) KST EP . X, Y , Z , X DOT , TOOT » ZDOT , I GOUT, JGOUT . KGOUT • PHI 

ISN 

0092 


URIT£(K,982)KE, JE.KB, JB,KA, JA.BE1, AL1,EV,P VOLTS 

ISN 

0093 

r 

STOP 

ISN 

009* 

3*0 

INT = 1 

ISN 

0095 


P H I SAV=PHI 

ISN 

0096 


SPE ED = 0. 

ISN 

0097 


PHI0LD=PH1 

ISN 

0098 

r 

IF(ME.GT.O) GO TO 350 

ISN 

0100 

L 

E = E E 

ISN 

0101 

r 

GO TO *00 

ISN 

0102 

u 

350 

CE=(A( JE) + FLOAT C2*KE“1“HE))/FL0AT(ME) 

ISN 

0103 


E=(1.+CE)/(1«-CE) 

ISN 

010* 


IFC XMACH.GT . 1. ) E=XKACH**2*(1.+CE)/(1. -CE) 

ISN 

0106 

r 

E=£ + AMAX1CPHI, 0.) 

ISN 

0107 

*00 

IFCE.LT.PHI) GO TO 1001 

ISN 

0109 

r 

SP£EO = SQRTCE'f*PHI) 

ISN 

0110 


XQOT=SPEED*SINA*CQS(BE) 

ISN 

Oil 1 


YDQT=SPEED*SINA*SINCBE) 

ISN 

0112 


ZDOT=SPEED*COSA 

ISN 

0113 


AL=ARCOS(COS A) 

ISN 

on* 


AL1=AL*180. /PI 

ISN 

0115 


BE1=BE*180./PI 

ISN 

0116 


EV=E*TEMP 

ISN 

0117 


P VOLTS =P HIS A V*T VOLTS 

ISN 

0118 


ZOLD=Z 

ISN 

0119 


KST EP = 0 

ISN 

0120 

r 

IFCNPRINT.NE.2.AND.NPRINT.NE.3) GO TO *90 



C PRINT INITIAL CONDITIONS OF TRAJECTORY 
r 

ISN 

0122 


WR I TECH, 982) KE , JE , KB, JB , KA, J A, B El , AL1 , E V , P VOLTS 

ISN 

0123 

982 

FORMAT C/1X.52HKE,JE, K8,JB, KA, JA, BETA, ALPH A , ENERGY* POTENTIAL 




1, /IX, 3(13, 12), 1PE22.8 ,*H 0£G,*X,E22.B ,*H DEG,8X,E16.*,2H V,*X, 



r 

2 £ 1 6. * , 2H V) 

ISN 

012* 

U 

WR ITECH,980) 

ISN 

0125 

980 

FORMAT ( 9X , 95HSTEPS X Y Z XDOT 



r 

1 YOOT ZDOT IG JG KG PHI) 

ISN 

0126 

L 

WRITEC M, 888 ) KSTEP ,X ,Y , Z ,XDOT , YD OT, ZDOT , IGOUT, JGOUT , KGOUT , PHI 

ISN 

0127 

888 

r 

FORMATC 9XjI5»lP6E11.3,3l6,E’11.3) 



L, 

c 

r 

TAKE A STEP 

ISN 

0128 

*90 

IF (KSTEP. EQ.O) GO TO 550 

ISN 

0130 

500 

CALL ORBIT 

ISN 

0131 


KSTEP=KSTEP + 1 

ISN 

0132 


IFCNPR INT.E0.3) WRITE(M,888) KST EP , X , Y , Z , X OOT ,YDOT , ZDOT , IGOUT , 




1 JGOUT, KGOUT, PHI 

ISN 

013* 


IF(KSTEP.LE.MSTEP) GO TO 550 

ISN 

0136 


WR ITE( M, 998) MSTEP 

ISN 

0137 

998 

FORMAT (/////IX, 9HM0RE THAN,I6,19H STEPS - HENCE STOP) 



> 


ISN 0138 
ISN 0139 
ISN 0141 
ISN 0143 


ISN 0145 
ISN 0147 


ISN 0149 
ISN 0150 
ISN 0152 
ISN 0153 

ISN 0155 
ISN 0156 

ISN 0157 
ISN 0159 
ISN 0160 
ISN 0161 


2 ISN 0163 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


0165 

0167 

0168 
0170 
0172 
0174 


0176 

0178 

0179 
0181 
0103 
0185 


0187 

0189 

0190 
0192 
0194 
0196 


C 

c 

550 


C 

560 


C 

570 


STOP 


IFCZ. EQ.O.. AND. ZDOT. LT.O. 

1 . AND . Y .GE. YPC13.AND. Y.LE.YPCJV33 GO TO 600 

IF (CX.LE.XX (1). AND.ZDOT.LT. 0-0. OR. 

1(X. GE.XXCII3. AND. ZDOT. LT.0.33GQ TO 600 

IFCCX.LE. XX (13. A NO. XDOT. LT.O.. AND. ZDOT. GT.O. 3. OR. 
1 (Y.LE.YY(l3.AN0.YD0T.LT.0.3.OR. 

2CX.GE. XXCII3.AND.XDOT .GE . 0 . . AND. ZDOT . GT . 0 . 3 . OR. 

3 CY.GE.YYC J J3 . AND. YDOT. GT.O. 3 . OR. 

4 (Z.GE.ZZCKK3. AND. ZDOI. GT.O. 33GO TO 700 


0 


538 


540 

C 


ISN 0198 


IFCSKPRFL.EQ.13 GO TO 538 

IF ( ( ( Y .LE.(YY(MBC3-.5*Z33.AND.(Y.GT •( Y YC HBC3-.5*ZZ(K.UK333 3 .OR. 
KCY.GE.C YY(MBD3+.5*Z33. ANO.CY.LT. C Y Y C M8D 3+ . 5*ZZ(KUK 3 33 3 
2.AND.X.GE.XX(ILX3.AN0.X.LE.XX(IUX33 GO TO 600 
CONTINUE 

IF (Z.NE.0..0R .ZOOT.GE.O.3 GO TO 540 
ZDOT =-ZDQT 

IF (NPRINT.EQ.33 WRITE(M,888) KSTEP , X , Y, Z , XDOT, YDOT, ZDOT, IGOUT , 

1 JGOUT .KGOUT ,PHI 
GO TO 590 
CONTINUE 

IF (KSTEP. EO. 03 GO TO 500 

PHIOLD=PHI 

CALL INTERP 

IFCIGOUT.LT.l.OR.IGDUT.GT.IIA.OR. JGOUT.LT.l.OR. JGOUT. GT.JJA.OR. 
1KGOUT. LT.l. OR. KGOUT. GT.KKA3 GO TO 330 
IFCNTIME.LT.1.0R.NTIME.GT.63 GO TO 580 

IFCNTIME.NE. 1 . AN D. NT I HE . N E . 2 3 GO TO 560 
X00TS = X00T**2 *■ PHIOLD-PHI 
IFCXD07S.EQ.0.3 XDOT=0. 

IFCXD0TS.GT.0..AND.X00T.NE.0.3' XDOT=SQRTC XD0TS3*SIGNC1 ., XDOT3 
IFCXDO TS. LT. 0. . AND. XDOT. NE. 0.3 XDOT* -XDOT 

IFCNPRINT.EQ.3.AND.XD0TS.LT.03 WRITECM,8883 KSTEP * X , Y, Z , XDOT , YDOT, 
1 ZDOT, IGOUT, JGQUT.KGOUT, PHI 

IFCNTIME.NE. 3. AND. NTIME.NE. 43 GO TO 570 
YD0TS=Y00r**2 + PHIOLD-PHI 
IFCYDOTS.EQ.0.3 YD0T=0. 

IF(YOOTS. GT.O.. AND. YOOT.NE. 0.3 YDOT=SQRT(YDOTS3*SIGNCl.,YDOT3 
IFC YOOTS. LT . 0. . AND. YDOT. NE. 0.3 YDOT=- YDOT 

IFCNPRINT.EO.3.ANO.YOQTS.LT.03 WRITECH.8883 KSTEP, X,Y,Z, XDOT, YDOT, 
1 ZDOT. IGOUT, JGOUT, KGOUT, PHI 

IFCNTIME.NE. 5. AND. NTIME.NE. 63 GO TO 590 
ZQQTS=ZD0T**2 + PHIOLD-PHI 
IFCZDOTS.EO.O.3 ZD0T=0. 

IFCZOOTS. GT.O. .AND. ZDOT. NE. 0.3 ZDOT=SQRT( ZD0TS3*SIGN(1 . ,ZD0T3 
IFCZDOTS.LT . 0. .AND. ZDOT. NE. 0.3 ZDOT=-ZDOT 

IF (NPRINT.EQ. 3. ANO.ZOQTS.LT .03 WRITECH, 8883 KSTEP,X,Y,Z, XDOT, YDOT, 
1 ZDOT, IGOUT, JGOUT, KGOUT, PHI 
GO TO 590 



ISN 

0199 

580 

ISN 

0200 

997 

ISN 

0201 


ISN 

0202 

887 

ISN 

0203 

C 

ISN 

020 4 

590 

ISN 

0205 


ISN 

0207 

1 

I SN 

0209 

C 

C PAR' 

ISN 

0210 

600 

ISN 

0211 


ISN 

0213 


ISN 

0214 


ISN 

0215 

C 

C PAR 
C 

ISN 

0216 

700 

ISN 

0217 


ISN 

0219 


ISN 

0221 


ISN 

0222 


ISN 

0223 


I SN 

0224 

720 

ISN 

0225 

740 

I SN 

0226 

C 

ISN 

0228 


ISN 

022 9 


ISN 

0230 


ISN 

0231 


ISN 

0233 


ISN 

0234 

C 

1 SN 

0235 

750 

ISN 

0237 


ISN 

0238 

889 

C 

ISN 

0239 

1002 

ISN 

0240 


ISN 

0242 


ISN 

0243 


ISN 

0244 


ISN 

0245 


ISN 

0246 


ISN 

024 7 


ISN 

0248 


ISN 

0249 

1000 

C 

C END 
C 

I SN 

0250 



FORMAT C/Z//1X,17HTR0UBLE - NTIME =,I3,19H = OUT OF RANGE 1-6} 
WRITECM, 887) KS T EP , X , Y , Z , XOOT , Y D OT , ZDOT , IGOUT, JGOUT , KGOUT , PHI , SAVE 
FOR MAT C 9X,I5,1P6E11.3,3I6,E11.3, , SAVE= , ,E18.10) 

STOP 

CALL INTERP 

IFC IGOUT.LI. l.OR.IGOUI.GT .IIA.OR.JGOUT.LT.l.OR.JGOUT.GT.JJA.OR. 
1KG0UT.LT. 1. OR. KGOUT. GT.KKA) GO TO 330 
IFCNPR INT.EQ.3) WRITECM, 888 ) KST E P , X , Y , Z , X DOT, YDGT , ZDOT, IGOUT, 

1 JGOUT, KGOUT , PHI 
GO TO 500 

ICLE IS ABSORBED' 

CONTINUE 

IFCNPRINT.NE.2.AND.NPRINT.NE.3) GO TO 1002 
FAT EC 1 )= ENO 1C 1) 

FATE(2)=EN01C2) 

GO TO 750 


CONTINUE 

IFCNPR INT.E0.1) GO TO 720 
IFCNPRINT.NE.2.AND.NPRINT.NE.3) GO TO 740 
FAIEC1)=END2C1) 

FATEC2)=EMD2C2) 

GO TO 740 

WRITECM, 98 2) KE , JE , K 8 , JB , KA , J A , B E 1 • AL1 , E V , P VOLTS 
NOESC=NOESC * 1 
IFCME.EQ.O) GO TO 750 

CSANGL=ZDOT/SQRTCXOOT**2+YDOT**2+ZDOT**2) 

XP0N=-2. *XMACH*SURTCE)*CSANGL - E - XKACH**2 
C0EFA=SPEEQ**2/FL0ATC NUMBER) 

IFCA8SCXP0N).GT.36.) GO TO 1000 
ADO =COEFA*EXPCXPON) 

DENS=DENS + AOO 

IFCNPRINT.NE.2.AND.NPRINT.NE. 3) GO TO 1002 

WRITECM, 889) FATE , KST E P , X , Y , Z , XOOT , YDOT » ZDOT , IGOUT » JGOUT, KGOUT ,PHI 
FOR MAT C IX, 2 A4, 15, 1P6E II. 3, 316, El 1.3) 

CONTINUE 

IFCMOSTPS.GE.KSTEP) GO TO 1000 

KE S-KE 

JES=JE 

KBS=KB 

JBS= JB 

KAS=KA 

JAS=J A 

MOSTPS=KSTEP 

CONTINUE 


FRACT=FLOATC NOE SC) /FLOAT C NUMBER) 



c 


ISN 

025 1 


ISN 

0252 

978 



C 

I SN 

0253 


ISN 

025 5 


ISN 

0257 

976 



C 

ISN 

0258 

800 

ISN 

0260 


ISN 

0261 


ISN 

0263 


ISN 

0264 

1001 

C 

ISN 

0265 

C 



c 



c 

ISN 

0267 


ISN 

0268 

972 

ISN 

0269 


ISN 

0270 

974 

T SN 

02 71 


ISN 

0272 



WRITECM.978) NOE SC , NU MB E R , FR ACT , E V, D E NS 

FORMAT C/1X,16HRATI0 ESCAPING = ,I5, 7H OUT 0F,I5,14H OR A FRACTION, 
1 F13.8,14H AT ENERGY E =,F13.8» 6H V QL T S , 4X . 6HC0ENS = , 1 PEI 4. 4, 1H ) ) 


IFCNPRINT.EQ.O) GO TO 800 
IF(ME.NE.O) WRITE(M, 976) 

FORMAT (IX, 80HDENS IS THE SUM OF AOD=SP EE D* *2*EXP (XP ON ) /NUMBE R OVER 
1 A HEMISPHERE OF DIRECTIONS//) 

I 

IFCME.EO.O) GO TO 1001 
C0EFE=2./C1. - C E ) ** 2 /FLO AT ( ME ) 

IFCXMACII.GT.l.) COEFE=COEFE*XMACH**2 

DENST=OENST + CQEFE*QENS 

CONTINUE 


IFCME.EO.O) DENST=SPEED**2*FRACT 


TRAJECTORY WITH MOST STEPS. PRINT K AND J INDICES. 


WRITE CM, 972) MOS T PS , KES , J E S ,K BS , JBS , K A S , JA S 
F0RMATC///1X, 15,3(13, 12), 29H =MOSTPS, KE.JE, KB, JB, 
WRI TECH, 974) XSAVE,YSAVE,PHISAV,DENST,PARTCL 
F0RMAT(/1X,26HAT niHENS IONLESS X.Y.PHI =, 3F12.6, 1H, , 
1 3 3H = NORMALIZED CURRENT DENSITY F0R.2A5//) 

RETURN 

END 


KA.JA) 

5X,1PE16.4, 



LtVfcL 21.8 C JUN 74 ) 


OS/360 FORTRAN H 


DATE 


ISN 0002 


COMPILER OPTIONS - NA ME = MA I N , OP T = 02 , L I NE C N T= 60 , ST ZE =0 000K, iinvnc _ 

SOURCE, FQCOIC.NOLIST, NODE CX, LOAD, MAP, NO EDIT, NOID, NO XREF 

SUBROUTINE 1NTERP 




ISN 0003 


ISN 0004 
ISN 0005 


ISN 0006 
ISN 0007 
ISN 0008 
ISN 0009 

ISN 0010 
ISN 0 0 L 1 
ISN 0012 


ISN 0013 
ISN 0015 
ISN 0017 

ISN 0019 
ISN 0020 
ISN 0021 
ISN' 0023 

ISN 0024 
ISN 0026 
ISN 0028 
ISN 0029 
ISN 0031 
ISN 0032 
ISN 0033 

ISN 0035 
ISN 0036 


ISN 0037 
ISN 0039 
ISN 0041 

ISN 0043 
ISN 0044 
ISN 0045 
ISN 0047 


C 

C 

C 


C 

c 

c 


10 

c 

100 

101 


102 

c 

103 

104 
C 

c 

c 

c 

c 

c 


INTERPOLATION WITHIN GRID 

COMMON/BK/IIM.IIP, J JM , J J P , KK , NT OT , I V . JV.II . J J.M, N , VP ( 3 0) , 

IX Y Z(20 80 ,39 , WC3Q.2Q, 101 , XP C 3 0) , XM( 10) ,YP(20),YMC10)»ZZC103t 
2XXC40),YY(30),ILX, I UX , KUK , MBC , MB 0 , V RF , NF P S , SKPRF L , SKP LST 
COMMON/O RO/XDOT , YDOT , ZDOT • X 1 » X 2 . Y1 , Y 2 • Zl , Z 2 , X , Y , Z , P HI , NT I ME , S A VE 
COMMON/ 1 NTER / I NT ,1 1 A , JJA ,KKA, I GOUT , JGOUT .KGOUT.XA, YA.ZA, 

1X1(30) ,YJ(20),ZK(10) 

IG0UT=0 
JGO UT = 0 
K. G 0 U T = 0 
NCH = 0 

X A = X 
Y A = Y 
ZA - Z 

LOCATE XA 

IFCXA.EQ.XICIIA)) IG=IIA-1 
IF(XA. EQ.XI ( 1 1 A ) ) GO TO 103 
IF(INT.NE.O) GO TO 100 


DO 10 1=2, IIA 
16=1-1 

IFCXA.LT. XICI)) 
CONTINUE 


GO TO 103 


IFCXA.GE .XICIG + D) GO TO 102 
IF(XA. GE.XICIG)) GO TO 104 
IG=IG-1 

IFCXA.LT. XICIG3) GO TO 101 

GO TO 103 

IG=IG*1 

IFCXA.GE. XICIGH33 GO TQ 102 

NCH = 1 
CONTINUE 

ACCEPT IF XI Cl G) LESS THAN OR EQUAL TO XA LESS THAN XICIG+1), 
LOCATE YA 

IFCYA.EQ.YJCJJA)) JG=JJA-1 
IFC YA. EQ.YJC JJA)) GO TO 203 
IFCINT . N E . 0 ) 


20 

C 


DO 20 J=2, JU A 
J G = J - 1 

IFCYA.LT. YJCJ3) 
CONTINUE 


GO TO 200 


GO TO 203 



! 


ISN 00*8 

ISN 0050 
ISN 0052 
ISN 0053 
ISN 0055 
ISN 0 0 5 6 
ISN 0057 

ISN 0059 
ISN 0080 


ISN 0061 
ISN 0063 
ISN 0065 

ISN 0067 
ISN 0068 
ISN 0069 
1 SN 0071 

ISN 0072 
ISN 007* 
ISN 0076 
ISN 0077 
ISN 0079 
ISN 0080 
. ISN 0081 
> 

5 ISN 0003 
ISN 008* 


ISN 0085 
ISN 0086 
ISN 0087 
ISN 0088 
ISN 0089 
ISN 0090 

ISN 0091 
ISN 0093 
ISN 0 0 9 * 
ISN 0095 

I SN 0096 
ISN 0090 
ISN 0099 
ISN 0 10 Q 


200 IFCYA.GE .YJC JG+1)) GO TO 202 
IFCYA.GE .YJC JG)) GO TO 20* 

201 JG = JG« 1 

IFC YA. LT .YJC JG)) GO TO 201 
GO TO 203 

202 J G = JG+ 1 

IFCYA.GE. YJCJG+l)) GO TO 202 
C 

203 NCH = 1 
20* CONTINUE 
C 

C ACCEPT IF YJUG) LESS THAN OR EQUAL TO YA LESS THAN YJCJ&n). 
C 

c 

C LOCATE Z A 

C 

IFCZA.EQ.ZKCKKA)) KG=KK A-l 
IFCZA.EQ.ZKCKKA)) GO TO 303 
IF(INT.NE.O) GO TO 300 

C 

DO 30 K=2»KKA 
KG=K-1 

IFCZA.LT.ZKCK)) GO TO 303 
30 CONTINUE 

C 

300 IFCZA.GE.ZKCKG + D) GO TO 302 
IFCZA.GE. ZKCKG)) GO TO 30* 

301 K & = KG- 1 

IFCZA.LT. ZKCKG)) GO TO 301 
GO TO 303 

302 KG=KG+1 

IFCZA.GE .ZKCK&+1)) GO TO 302 
C 

303 NC II = 1 
30* CONTINUE 
C 

C ACCEPT IF ZKCKG) LESS THAN OR EQUAL TO ZA LESS THAN ZKCKG*1). 

C 

C LOCATE LINE AND BOX 

C 

XI^XIC 1G) 

Y1=YJC JG) 

Zl = ZKCKG) 

X2 = XIC 1G+1) 

Y2=YJC JG + 1) 

Z2-ZKC KG + 1) 

C 

IFCX.NE.X1.QR.XDOT.6E.O.) GO TQ *00 
IG-TG-1 
X 2 = X 1 
X1=XICI&) 

C 

*00 IFCY.NE.Y1.0R.Y00T.&E.0. ) GO TO 500 

JG= JG-1 
Y2= Y 1 
Y1=YJ( JG) 

C 

500 IFCZ.NE.Zl.OR.ZDOT.GE.O. ) GO TO 600 



ISN 0101 


1SN 0103 
TSN 0104 
I SN OIOS 

ISN 0106 
1 SN 0107 
ISN 0108 
TSN 0109 
ISN 0110 
ISN 0111 


c 

600 


KG=KG-1 
Z2 = Z1 
Z1 = ZKXKG) 


PHI=VV(IG»JG»KG) 

IGUUT = 1G 

JGOUT= JG 

KGOUT=KG 

RETURN 

END 


COMPILER OPTIONS - NAME = M A I N , 0 P T = 02 • L I NE C NT = 60 , S 1 1 E = 0 000K , 


ISN 0002 


ISN 0003 


ISN 0004 


1 SN 
ISN 
ISN 
ISN 
T SN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 

ISN 

ISN 


0005 

0006 
0007 
0009 
'0010 
0012 
0014 
0016 
0018 

0019 

0020 

0021 

0022 


ISN 0024 


I SN 
ISN 
> ISN 
£ ISN 
ISN 


0025 

0026 
002 7 
0028 
002 9 


C 

C 

c 


997 


990 

991 

992 

993 
C 


NR M fc = Mfl 1 N, UK I = Ud • L I NCL N I -OU f 3 1 tn-UUUUJ • 

SOURCE, EBCDIC, NOLIST,NOQECK»LQAD,MAP»NOEQIT»NOID*NQXREF 
SUBROUTINE POWER 

CURRENT DENSITIES AND POWER LOSS 

COMMON/BK/IIM.ITP, JJM, JJP,KK,NTOT,IV, JV,II,JJ,H,N,VP(30}, 
lXYZ(2080»3)»VVC30»20,10),XP(30)tXM(lO),YP(20),YM(lQ)»ZZ(lO!» 
2XXC403 ,YYC303 . I LX . IUX , KUK , MBC, MB D, VRF , NFP S , SKPRFL , SKPLST 
CUMMON/CP/NPRINT.NPrS.MA.MB.MEf KMAX.XPT.YPT.ALl.BEl.EV.SHACH, 

1 TVQLTS.CUR, XMETER 

CO MMON/CO/P VOLTS , X M A CH . 0 E NS T , NN , P AR TCL ( 2 3 . P ART1C 23 . P ART2C 23 
DIMENSION A C 2 ) 

IFCNPTS.EQ.O.OR.MA.EQ.O! WRITECM.997) X P I . YPT , AL1 . BE 1 » EV 
FORHATC/IX, 9HX ANO Y = , 2F10. 5, 20X, 19 HALP HA, BETA, ENERGY =,3F20.5) 
IFCNPR INT.EQ.03 WRir£(M,990! 

IFCNPRINT.E0.1) WRI T EC M, 9913 
IFCNPRINT.E0.23 U/R I T E C M , 9 9 2 3 
1FCNPRINT.E0.3) WRITECM.993) 

FORMAT C/38H NPRINT=0 MEANS NO TRAJECTORY PRINTINGS 
FORMAT C/53H NPRINT = 1 PRINT INDICES OF ESCAPING TRAJECTORIES ONLY! 
F0RMATC/56H NPRINT=2 PRINT FIRST ANO LAST STEPS OF ALL TRAJECTORIE 
ISO 

FORMAT C/52H NPRINT = 3 MEANS PRINT EVERY STEP OF ALL TRAJECTORIES! 


C 

C 

C 


IFCTVOLTS.EQ.O.) RETURN 
X MACHISM ACM 


200 
C 


NON-DIMENSIONALIZE THE POTENTIAL DISTRIBUTION. THEN RESTORE AT END. 

DO 200 K = 1 . K K. 

DO 200 J=1,JJ 

00 200 1 = 1,11 

VV(I,J,K 9=VVCI, J.IO/TVOLTS 
CONTINUE 


C 

C 

C 

C 


DEFINE THE PANEL POINTS AT WHICH THE CURRENT AND POWER IS EVALUATED 
CASE OF A SINGLE POINT 


ISN 

0030 

C 

C 

f 

IFCNPTS.EQ.O.OR.MA.EQ 
CASE OF MULTIPLE POINTS F 

ISN 

0032 

L. 

JV M = 1 

ISN 

0033 


I V M = 1 

ISN 

0034 


IF(JV.GT.l) JVM=JV-1 

ISN 

0036 


IFCIV.GT.I) IVM=IV-1 

ISN 

0039 


N A= 0 

ISN 

00 3 9 


N AR E AS = I VM* JVM 

ISN 

004 0 


TPOWER=0. 

ISN 

004 1 


TCURNT = 0 . 

ISN 

0042 


T ARE A = 0. 

ISN 

0043 


.NN = 0 

ISN 

004 4 


DO 500 J = 1 » J V M 

r sn 

004 5 


DO 500 1=1, IVM 

ISN 

0046 


N A=NA+ 1 


- • ■ 

date 


r sn 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

rSN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


00A7 
00 A B 

0050 

0051 

0052 

0053 
005 A 

0055 

0056 

0057 
005 3 

ms 
0061 
0062 
0063 
00 6 A 

0065 

0066 

0067 

0068 
0069 
007 1 

0072 

0073 
00 7 A 

0075 

0076 

0077 


> 

to 

o 


250 


260 

270 


IFCNPTS.EQ.O.OR.MA.EQ.O) GO TO 250 
PO = 0 
CU. = 0 

AC1)=-1./SQRTC3.) 

AC2)=-AC1) 

GO TO 260 
CONTINUE 
JXMAX=1 
J Y M AX = 1 
K M A X= 1 
GO TO 270 
JX HAX = 2 
JY M AX = 2 
CONTINUE 

00 AO 0 K Y = l» KMAX . 

00 AO 0 KX = 1.KHAX 
00 A00 JY = 1 , J YM AX 
00 A00 JX= 1 » JXMAX 
NP = NP 4-1 
NN-=NN + 1 

IFCNPTS.EQ.O.OR.MA.EQ.O) GO TO 300 
CX=CACJX) * FLOAT C2*KX - 1 - K M A X ) ) /FLO AT CKMAX) 
CY=CACJY) + FLOAT C2*K Y - 1 - KM AX))/FL0AT C KMAX) 
XPT=CXPCI+1)-XPCI))/2.*CX > (XPCI+l)+XP(I))/2. 
YPT = CYPC J+1)-YPCJ))/2.*CY + C YP( JUUYPC J))/2. 

coef = cxpch-i)-xpci))*cypcj+i)-ypcj>) 

AREA = COEF * XHETER**2 
COEFM = AREA/A. /FLOAT CKMAX **2) 






\ 


. T » 

/,rx«r - Nt! 


ISN 

0078 

300 

ISN 

0079 


ISN 

0080 


ISN 

0081 


ISN 

0083 


ISN 

008 A 


ISN 

0085 


ISN 

0086 


ISN 

0087 


ISN 

0088 

C 

ISN 

00 89 

995 



C 

ISN 

0090 


ISN 

0092 

9 9 A 



C 

ISN 

0093 


ISN 

00 9 A 


ISN 

0095 

988 


MULTIPLY BY LOCAL POTENTIAL TO 


C COMPUTE EACH CURRENT DENSITY AND 
C EVALUATE POWER DENSITY 
CALL DEN 

DENCUR=OENST*CUR 
POUDEN=PVOLTS*OENCUR 
IFCMA.EQ.O) GO TO 600 
XPTM=XPT*XMETER 
YPTM=YPT*KMETER 
XPM=XPCI)*XM£TER 
XPPM=XP(IH)*XMETER 
YPM = YP CJ)*XMETER 
YPPM=YP( J+1)*XMETER 

FORMAT C6X.12HAT POINT NO..I3.10H, WITH X .=, F10. 5. 13H METERS . Y = 
ltF10.5 « 27H METERS, ANO COEFFICIENT =,F10.5,1AH SQUARE METERS) 

IFC'NPTS.GT.O.ANO.MA.GT .0) WRITECM, 99A) N A . XPM.XPPH, YPM , YP PM 
FORMAT C /5X.16H IN SUB-AREA NO. , 13, IX , 17HBOUNDED BY X IN C, 

1 F10.5.3H T0.F10.5, 9H) METERS ,, AX , 1 3HAND BY Y IN C, 

2 Fl’0.5 ,3H TO,F10.5, 8H) METERS) 

WRITE CM, 995 ) NP , X PTM , YP TH , COEFM 
WRIT EC M, 988) PVOLTS,DENCUR,POWDEN,PARTCL 

FORMAT C6X,53HTHE VOLTAGE, CURRENT DENSITY, AND POWER DENSITY ARE = 
1/6X * IP £ .1 6. A , 6H VOLTS, AX, £16. A, 23H AM P/C SO- METER) , AND , E16.A, 

2 2 A H WATT/CSQ-METER), F0R.2A5//) 

IFCNPTS.EQ.O ) GO TO 600 


■t L 
3^ * 


L TKS C 1 - IT 


ISN 0096 



ISN 

009B 


ISN 

0099 


ISN 

0100 

400 

ISN 

0101 


ISN 

0102 

C 

I SN 

0103 


I SN 

0104 


ISN 

0105 

9 B 4 

ISN 

0106 

986 



C 

ISN 

0107 


I SN 

0108 


ISN 

0109 


ISN 

0110 

500 



C 

ISN 

0111 


ISN 

0112 

982 

I SN 

0113 


T SN 

0114 


ISN 

0115 


ISN 

0116 

980 


c 

c 

c 

c 


CU = CU + COEFM*DENCUR 

PO = PQ + COEFM * PQWOEN 

CONTINUE 

AVCD = CU/AREA 

A V P 0 = PO/AREA I 


WRITEC M, 986) N A , CU , P 0 , P AR TCL 
WRITECM, 984) NA,AREA,AVCD,AVPD 

FORMATC IX, IB H IN SUB-AREA NUMBER, I3*BH OF AREA, 1PE16.4, 15H SQUARE 

1 METERS, /52H THE AVERAGE CURRENT DENSITY AND POWER DENSITY ARE = , 

2 E 16.4, 19H AMP/CSQ-METER ) AND , E 1 6. 4 , 1 6 U WATT/CSQ-METER)) 

FORMATC /1X.18HIN SUB-AREA NU MB E R , 1 B , 2 8H THE CURRENT AND POWER ARE 
1 =,1PE16.4»12H AMP, AND.E16.4 ,HH WATTS, FQR.2A5) 


T AREA=TARE A + AREA 
TCURNT = TCURNT > CU 
TPOWER = TPOWER + PO 
CONTINUE 


WRITECM, 982 ) TCU R N T , T P 0 WE R , P A RT C L 

FORMAT C///1X,34HT0TAL CURRENT AND POWER LOSS ARE =,1PE16.4, 

1 1 2 H AMP, AN0.E16. A, 13H WATT, F0R,2A5) 

AVCD=TCURNT/TAREA 

AVPD=TPOW£R/TAREA 

WRITECM, 980) T AR E A , A VCO , A VP D 

FORMAT C/IK. 26HWITH A TOTAL PANEL AREA OF , 1 PE 16. A , 1 5H SQUARE METERS 
1, /IX, 51 H THE AVERAGE CURRENT DENSITY AND POWER DENSITY ARE = , 

2 E16.4.19H AMP/CSQ-METER) AND , E 1 6.4 , 1 6H WATT/CSQ-METER)) 


RESTORE POTENTIAL DISTRIBUTION TO DIMENSIONAL VALUES 



ISN 

0117 

600 

CONTINUE 

> 

ISN 

Oil 8 


DO 700 K = 1 , X K 

DO 

ISN 

0119 


DO 700 J=1,JJ 


ISN 

0120 


DO 700 1=1,11 


I SN 

012 1 


VV(I,J,K)=VV(I»J,K)*TVOLTS 


ISN 

0122 

700 

CONTINUE 


ISN 

0 3 2 3 


RETURN 


ISN 

0124 


END 


A22 


COMPILER OPTIONS - NAME = M A I N, OP T =02 , LINE C NT= 60 , S1ZE = Q OOOK , 

SOUR C E, E BC 01 C, NOLIS T , NO D ECK, LOAD .MAP, NOEDIT, NO ID, NOXREF 
ISN 0002 SUBROUTINE LISICLST.IP) 

ISN 0003 C0MM0N/BK/IlM,IIP f JJM,JJP,KKtNT0T,IV,JV,II.JJ,M,N,VPC3 0'>. 

lXYZC20aO,3),VVC30,20,10),XPC30),XMClO),YPC20),YHC10),ZZC10>» 
2XXC40) , Y YC30) , I L X , I UX , KUK , M8C , MB D , VR F , NF P S , SKPRF L, SKPLST 
ISN 0004 COHMON/FLD/X C2080, 2) .COEFC2080, 7) .INDXC2080, 6), SKPCO 

ISN 0005 DIMENSION KQU T C 5 ) , XOU T C 5 ) , YOU T C 5 ) , ZQUT C 5 ) 

TSN 0006 DO 500 LINE=1.60 

ISN 0007 00 200 NP=1, 5 

TSN 0008 KP=LINE * CNP-1)*60 + CIP-1) * 300 

ISN 0000 IFCKP.GT.NTOT.AND.NP.EQ.il RETURN 

TSN 0011 IFCKP .GT. NTOT) GO TO 300 

ISN 0013 N M A X = N P 

ISN 0014 KOUTCNP) = KP 

ISN 0015 IFCLST.EQ.l) XOUTCNP) = XCKP.l) 

ISN 0017 IFCLST .EQ. 2) XOUTCNP) =XYZCKP,1) 

ISN 0019 IFCLST .EQ. 2) Y OUTC N P )=X YZ CKP , 2 ) 

ISN 0021 IFCLST .EQ. 2) ZOUTCNP) =XYZCKP,3) 

TSN 0023 200 CONTINUE 

ISN 0024 300 GO TO C 400 , 4 50 ) , LST 

ISN 0025 400 WRITECM.1000) C K OUT C NP ), X OUTC NP ) , NP=1,NMAX) 

ISN 0026 1000 FORMAT C5CI8»F16.8)) 

ISN 0027 GO TO 500 

Hi! §§?§ 3000 ^;| C T ?^?g;> F «»^^' > >-X^tN P 3,,0UrCNP),Z0UT CN |.),NP = l,NKAX) 

TSN 0030 500 CONTINUE 

ISN 0031 RETURN 

ISN 0032 END 



US/ 30U rUKIKAN H 

COMPILER OPTIONS - NAME = M A I N , 0 P T =02 , L I NE C NT= 6 0 , S I ZE =0 000K , 

SOURCE, EBCDIC, NOLIST.NQDECK, LOAD, MAP, NOEOIT.NOID.NQXREF 
; SUBROUTINE RELAX 

C POINT-SUCCESSIVE 0 VE R R E L AX AT I ON HETHOD 

COMM 1)N/BK/IIM,1IP»JJM»JJP,KK,N TOT «IV»JV,II,JJ,M»N»VP(30)» 

1XYZ (2080,3), V V ( 3 0 , 20, 10), XP(30),XM(10) «YP(20)»YM(10)»ZZ(10)» 
2XX(A0),YY(30),ILX, I UX.K UK, MBC, MB D.VRF, NFPS.SKPRFL.SK PLST 
: COMMON /FLO/X C208 0, 2) , CO E F C 2 0 8 0 , 7 ) , I NO X C 2 0 8 0 , 6) , S KPC 0 

0 ME GA= 1.9 

, EPS = l.E-3 

ITMAX=2000 
I T R = 0 

i JPROLD=0 

I G 0 = 1 

200 ITR=ITR+1 
DELTA M=0. 

i DO 500 N = 1 , NTQT 

IF(X(N,2).EQ.1)G0 TO 500 
, X1=X(N,1) 

C 

r FN=C0EF(N,1)/C0EF(N,7) 

] FS=C0EF(N,2)/C0EF(N,7) 

I FE=C0EF(N,3)/COEF(N,7) 

) F'J = COEF(N,4)/COEF(N, 7) 

FU=COEF(N, 5)/C0EF(N, 7) 

! F0=C0EFCN,6)/C0EF(N,7) 

c 

l NN=INDXCN,1) 

, NS=INQX( N, 2) 

; NE=INDX(N,3) 

» NW=INOX( N, 4 ) 

NU=IN0X(N,5) 

) no=inox<:n,6) 


SUM=0. 

IF(NN.GT.O) SUM 
IF(NS.GT.O) SUM 
IFCNE.GT.O) SUM 
IFCNW.GT.O) SUM 
IF(NO.GT.O) SUM 
IF(NU.GT.O) SUM 


SUM+FN*XC NN, 1) 
S'JM«-FS*XCNS,1) 
SUM+FE*XCNE,1) 
SUM«-FW*X(NW, 1 ) 
Sl)M + FO*X(ND, 1) 
SUM*-FU*X(NU,1) 


X (N , 1 ) = OMEGA*SUM+(l.-QMEGA)*Xl 
DELTA = ABS(X(N,1)-X1) 

IFCm(Xn.&T.l.E-10) 0ELTA = ABSC(XCN,1)-X1)2X1) 

IFCDELTA .GT. OELTAM) OELTAM=DELTA 

CONTINUE 

IF(ITR.GT.ITMAX) WRITEC M, 8888) ITR 
IF( ITR. GT. UMAX) GU TO 700 

FORMAT (////I OH MORE THAN, I4.11H ITERATIONS) 
IPR=ITR/500 

IF(IPR.LE.IPRQLO) GO TO 600 

IPROLD=IPR 

GO TO 800 

IFCOELTAM.GT.EPS) GO TO 200 


ITERATION FINISHEO. PRINT AND EXIT. 


A 24 


C 

I S N 006 1 700 I GO-2 

ISN 0062 BOO NFP P = (NTOT/300) +1 

ISN 0 0 6 3 LH) 900 IP=1»NFPP 

ISN 0 0 64 WRITE CM, 777 7 0 ITR, EPS, DELTAM, OMEGA 

ISN 0065 7777 F URM A T ( 1 5111 S GLUT ION AFT E R, 16 , 2X , 25HI TE R AT I ONS WITH TOLERANCE, 

1 Fl2.H»0X« 18 H MAXIMUM 0 1 F F E R E N C E , F 1 2 . 8 , B X , 6 HQ MEG A= , F 8 . 5 ) 

ISN 0066 CALL LISTCl.IP) 

TSN 0067 900 CONTINUE 

ISN 0060 U 6-0 — TO— <-6-0 0-.-1-0 0 OIvTC 0 (F C I 0 • £ ? • 1 ) Of f 0 L O D 

ISN -00 69 1 0 00 RETURN 

ISN 0070 END 



A25 


LEVEL 21.8 ( JUN 74 ) 


□S/360 FORTRAN H 


TSN 
I SN 


I SN 
! SN 
1 SN 
1 SN 
1 SN 
1 SN 
I SN 
ISN 
I SN 


COMPILER OPTIONS - NAME = M A I N , 0 P T = 0 2 , LI N E C N T = 60 , S I ZE = 0 GOOK, 

source,eocoic,nolist,nqdecic,load,map,noeoit,noid,noxref 

0002 SUBROUTINE FINO(I,J,IO 

00 0 3 COMMON/BK/ II M , 1 1 P , J JM , J JP , KK , NTOT . I V . J V , 1 1 . J J, H , N , VP ( 3 0) , 

lXVZ(2080*3)*VVC30»20,10)*XP(30)»XM(10).YP(20)«YMClO)jZZ(lQ)» 
2XXC40) » Y YC 30 ) . I L X , I UX , KUK , MBC , MB 0 , V RF , NF P S , SKPRFL , SUP LS T 

0004 IT J J=T l* JJ 

0005 K=N/IIJJ+1 

0006 IFCK .GE. 2 .AND. MOO(N,IIJJ) .EO. 0) K = K-1 

0008 N K I J- N - IIJJ*(K-1) 

0000' J=NKIJ/IIH 

00 10 I F ( J .GE. 2 .ANQ. MQDCNKIJ ,11) .EQ. 0) J = J-1 

0012 I = N K I J - II*(J-1) 

0013 RETURN 

0014 EMU 



LLVtL 21.8 ( JUN 74 ) 


QS/360 FORTRAN H 


DATE 


> 

to 

OT 


TSN 0002 
1SN 0003 


TSN 0004 
ISN 0005 


COMPILER OPTIONS - NAME = M A I N , 0 P T = 0 2 . L I NE C NT = 60 , S T 2 E = 0 000K . 

SOURCE, EBCDIC, NOLI ST, NODECK , LO AD , M A P ,NOEDI T , NO ID , NO XREF 
SUBROUTINE ARRAY 

COMMON/ BK/IIM, IIP, JJM,JJP,KK,NTOT ,IV,JV,II,JJ,M,N»VP(30)» 
1XYZC20U0 ,3) , V V C 30 ,20 ,10) »XPC30) ,XM( 1 6 ) ,YP( 20),YMC10) ,ZZC10)» 
2XXC40),YY(30),ILX»IUX»KUK,MBC»M&0,VRF»NFPS,SKPRFL,SKPLST 
COMMON/FLD/XC2080»2),COEFC2Q80,7),INDXC2080»6)«SKPCO 
COMMON/CCC/CN,CS,CE,Cy,CU,CO,CC,NN,NS,NE,NW,NU,ND 


TSN 0006 
ISN 0007 
ISN 0008 
ISN 0009 
ISN 0010 
ISN 0011 
ISN 0012 

ISN 0013 
TSN 0014 
ISN 0015 
ISN 0016 
ISN 0017 
ISN 0018 
ISN 0019 
ISN 0021 
ISN 0022 
ISN 0023 
ISN 0024 
ISN 0025 


C COEFFICIENT ARRAY = COEFCN, 7), WHERE 

C C0EF(N,n = CN CN0RTH=*Y NEIGHBOR! 

C COEFCN, 2)=CS (SOUTH=-Y NEIGHBOR) 

C C0EF(N,3)=CE CEAST=+X NEIGHBOR) 

C COEFCN, 4 ) = C W C WEST=-X NEIGHBOR) 

C COEFCN, b) = CU C UP = + Z NEIGHBOR) 

C COEFCN, 6)=CD C DOWN=-Z NEIGHBOR) 

C COEFCN, 7)=CC C = CENTRAL POINT) 

C 

C SAVE COEFFICIENTS AND INDICES 

C 

COEFCN ,1)=CN 
COEFCN , 2 ) = C S 
COEFCN, 3)=CE 
COEFCN ,4 ) = CW 
COEFCN, 6) = CU 
COEFCN, 6) = CD 
COEFCN, 7)=CC 
C 

I ND X C N , 1 ) = NN 
INOXCN, 2) = NS 
I N 0 X C N , 3 ) = N E 
INOXCN , 4 ) = NW 
INOXCN , 5) = NU 
INOXCN, 6)=ND 
IFCSKPCO.EQ. 1) GO TO 20 

WRITECM, 1000) NO,CO,NS,CS,NW,CW,N,CC,NE,CE,NN»CN,NLI,CU 
1000 F0RMATC/7C1X , 1HC , 1 4, 2H) = , 1 PEI 0 . 4) ) 

20 CONTINUE 
RETURN 
ENO 



A 27 


ISN 0002 
ISN 0003 


ISN 

ISN 

ISN 

1 SN 

I SN 

ISN 

T SN 

ISN 

I SN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

I SN 

TSN 

ISN 

ISN 

ISN 

I SN 

ISN 

ISN 

TSN 

ISN 

ISN 

ISN 

ISN 

ISN 

I SN 

T SN 

TSN 

T SN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


0004 

0005 

0006 
0007 
000 B 
0009 
0011 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
0025 

0027 

0028 

002 9 

0030 

0031 

0032 
00 33 
00 3 4 

0035 

0036 

003 7 
0039 

004 1 
0042 
004 3 

0045 

0046 

004 7 
00 4 8 

0049 

0050 
0052 

005 3 
0054 
00 55 
0056 
005 7 


QS/360 FORTRAN H 

COMPILER OPTIONS - NAME= M A I N, OP T = 02 , LI NE C NT = 60 , S 1 1 E = 0000K, 

SUBROUTINE cSDCMP^cia5 NaLIST ’ r<OC ’ 6CiC * LOAD ' HAP ’■ NOeOIT,!JOIO ’ MOXREF 

2XX(40)*YY(3&.) J lLXf IUXiKUK#MBC»M8DfVRF|NFPS»SK?RFL«S!CPLST f 

NF^ , if^^5 N,c h^;S^,5^ : ^- c 9!l!"- Ns:N E'»"-NuVND’ SICPLST 

NFC IX, J X ,KX1 - IX+II*( JX-13 + II*JJ*(KX-13 

A - U • 

C = 0. 

CALL FINDCI.J.K3 
IFCI .EQ. 1) GO TO 100 
IF C I .EQ. II 3 GO TO 200 
NH = NF( I+lt J» K) ■ 


XYZ(NL.l) 


- XYZ(N,13 


NL=NF( I-l.J.K) 
DX=XYZ(NH,1) - 
GO TO 300 
100 NH=NF(2, J » K 3 
DX=XYZ(NH, 1) 

GO TO 300 
200 NL=NF(1I-1, J.K) 

0X=XYZ(N,1) - XYZ(NL.l) 
300 CONTINUE 

IFCJ .EO. 13 GO TO 400 
IFCJ .EQ. JJ3 GO TO 500 
NH=NF(I,J+l,K3 
NL=NFC1, J-l.K) 
DY=XYZCNH,2) - XYZCNL.23 
GO TO 600 

400 NH=NF(I,2,K3 

DY = XYZ (NH* 23 - XYZCN.23 
GO TO 600 

500 NL=NF( I , JJ-l.K) 
n _ 0Y = XYZ(N,2) - XYZCNL.23 

60 0 A = DX*D Y/4. 

1F(MP .EQ. n GO TO 700 
IFCMP .EO. 23 GO . TO 800 
RETURN 
700 N 0= 0 

IFCK .EO. KK 3 RETURN 
NH=NF( I, J.K+13 
N U = N H 

OZ=XYZCNH,33 - XYZCN.33 
GO TO 900 
.800 N0=0 

IFCK .EQ. 13 RETURN 
NL-NFC I, J.K-I3 
ND=NL 

boo UDIY”-” • XY2(NL ’ 3 > 

RETURN 

END 


UATE 


LEVtL 21.6 ( JUN 7 4 ) 
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DATE 


> 

to 

00 


COMPILER 


1 SN 

0002 


I SN 

000 3 


I SN 

0004 


I SN 

000 5 


I SN 

0006 


I SN 

000 7 


1 SN 

00 OH 


T SN 

0009 


] SN 

00 11 


T SN 

0013 


1 SN 

0014 


1 SN 

0015 


I SN 

0016 


1 SN 

001 7 

100 

T SN 

00 IB 


1 SN 

00 19 


1 SN 

0 0 2 0 

200 

! S N 

00 2 1 


T S N 

002 2 

300 

1 SN 

00 2 3 


I SN 

002 5 


1 SN 

002 7 


] SN 

00 2 8 


1 SN 

0029 


1SN 

0030 


ISN 

003 1 

400 

I SN 

00 32 

T SN 

0033 


ISN 

0034 

500 

ISN 

0035 


T SN 

0036 

600 

ISN 

0037 

I SN 

00 3 9 


ISN 

0 0 4 1 


ISN 

0042 

700 

ISN 

0043 

I SN 

00 4 5 


TSN 

004 6 


ISN 

004 7 


TSN 

0048 


1 SN 

0049 

800 

TSN 

005 0 


ISN 

0052 


ISN 

005 3 


ISN 

00 5 4 


TSN 

005 5 

900 

! SN 

005 6 


1 SN 

00 5 7 



OPTIONS - NAME= M A I N , 0 P T = 0 2 . L I N E C N T = 6 0 , S 1 Z E = 0 OOOK , 

SOURCE .EBCDIC, NOLIST, NOD ECK,LO AD, MAP,NOEDIT, NO ID. NOXREF 
SUBROUTINE CNSCMP,C,AI 

COMMON/fiK/IIM, I IP , JJM , J JP , KK , N TOT , I V, JV.II , J J, M , N , VP C 3 OT , 
1XYZ(2080»3),VV(30»20,10},XP(30),XM(10),YP(20),YM(10),Z7(10), 
2XXC 40) , Y Y(30 ) , ILX, IUX, KUK , MBC,MBD, VRF, NFPS , SKPRFL.SKPLST 
COMNON/CCC/CN,CS,CE,CW,CU,CD,CC f NN,NS,NE,NW,NU,NO 
NF(IX,JX,KX) = IX*II*(JX -1)4-11 *JJ*CKX-1) 

A =0 . 

C- 0 . 

CALL F INDC I, J, K) 

IFCI. E0.1) GO TO 100 
IF( T.E (J. II) GO TO 200 
NU-NFC 1+ 1, J, m 

N L = N F ( I-lfJflC) 

DX = XYZCNI!,n - XY2C.NL, 1) 

GO TO 300 
NU=NF(2, J,K) 

0 X = X Y 7. ( N H * 1) - XYZ(N.l) 

GO TO 300 
NL=NFCII«1,J,K3 
0X=XYZ(N,1) - XYZCNL.l) 

CONTINUE 

IFCK.EQ.l) GO TO 400 
IFCK.EQ.KIO, GO TO 500 
NU = NIFC f, J.K4-1) 

NL=NF(I, J,K-1) 

0Z=XYZ(NH,3) - XYZCNL.3) 

GO TO 600 
NH=NFCI, J, 2) 
l)Z=XYZ(NH,3) - XYZCN.3) 

GO TO 600 
NL=NF(I, J.KK-1) 

DZ = XYZ(N ,3) - XYZCNL.3) 

A=0X*DZ/4. 

IFCMP.EO.l) GO TO 700 
IFCMP.EQ.23 GO TO 800 
RETURN 
NN = 0 

IFCJ. EQ.JJ) RETURN 
NU = NFC I, J*l, K) 

NN = NH 

0Y=XYZ(NH,2) - XYZCN.2) 

GO TO 900 
NS = 0 

IF(J.EU.l) RETURN 
NL=NF(I, 

N $ =NL 

DY=XYZ(N,2) - XYZ(NL,2) 

C = A /D Y 
RETURN 
ENO 
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OS/ 360 FORTRAN H 

COMPILER OPTIUNS - NAME = M A I N, OP T =0 2, LINE CNT = 60 . S T ZE - QGOOK , 

SOURCE, EBCDIC, NOLIST, NOOECK, LOAD, MAP, NOEOIT, NO ID »NOXREF 

ISN 0002 SUBROUTINE CEUCMP ,C, A) ■ 

ISN 0003 COMMON/BK/IIM, IIP, JJM, J J P , KK , NT 0 T , I V , J V , I I , J J,M ,N , VP C3 0) , 

1XYZC2080 ,3), VVC30.20 , 10) ,XPC3Q3, XMC103 , Y P C 20 3 , YM C 103 , Z ZC 1 0 3 , 
2XXC40).YYC30),ILX,IUX,KUK,MBC,M8D,I/RF,NFPS,SKPRFL f SKPLST 
ISN 0004 COMMON/CCC/CN,CS,CE,CW,CU,CD,CC,NN,NS,NE,NW,NU,ND 

ISN 0005 NFC IX, JX ,K.X3 = IX+ Il*CJX-13 + 1 1* j j* c k.x - 1 3 

ISN 0006 A=0. 

ISN 0007 C=0. 

ISN 0008 CALL FINDCI.J.K3 

ISN 0009 IF(J.EQ.l) GO TO 100 

ISN 0011 IFCJ .EQ. JJ3 GO TO 200 

ISN 0013 NH = NF(I,JH,K3 

ISN 0014 NL=NF(I, J-l.K) 

ISN 0015 QY=XYZCNH» 2) - XYZCNL.2) 

ISN 0016 GO TO 300 

ISN 0017 100 NH=MF(Ij2,K) 

ISN 0018 DY=XYZCNH,2) - XYZCN.23 

ISN 0019 GO TO 300 

ISN 0020 200 NL=NFC1. JJ-1.K3 

ISN 002 1 DY=XYZCN,2) - XYZ(NL,2) 

ISN 0022 300 CONTINUE 

ISN 0023 IFCK .EQ. 13 GO TO 400 

ISN 0025 IFCK .EQ. KK) GO TO 500 

1 SN 0027 NH=NFC I, J, K.+ 13 

ISN 0028 NL=NFCI, J,K-13 

ISN 0029 DZ=XYZCNH,3) - XYZCNL.33 

ISN 0030 GO TO 600 

ISN 0031 400 NH=NFCI,J,2) 

ISN 0032 DZ=XYZCNH,33 - XYZ(N,33 

ISN 0033 GO TO 600 

ISN 0034 500 NL=NF(I, J.KK-13 

ISN 0035 DZ=XYZCN,33 - XYZCNL,3) 

ISN 0036 600 A=DY*DZ/4. 

ISN 0037 IFCMP .EQ. 13 GO TO 700 

ISN 0039 IFCMP -EQ. 23 GO TO 800 

ISN 0041 RETURN 

ISN 00 42 700 NE = 0 

ISN 0043 IFCt.EQ.II3 RETURN 

ISN 0045 NH=NFC 1+ 1 , J, K ) 

ISN 0046 N E = N H 

ISN 0047 DX=XYZCNH, 13 - XYZCN.l) 

ISN 0048 GO TO 900 

ISN 0049 800 NW=0 

ISN 0050 IFCI .EQ. 13 RETURN 

ISN 0052 NL=NFCI-1.J.K) 

ISN 0053 NW = NL 

ISN 0054 DX = X YZ (N ,13 - XYZCNL.l) 

ISN 0055 900 C= A/DX 

ISN 0056 RETURN 

ISN 0057 ENO 


DATE 


UJ/ JUU 


rui\UAN n 


> 

to 

© 


ISN 0002 


COMPILER OPTIONS - NA ME = MAIN,DPT = 02»LINECNT = 60f51/. E-0000K, 

SQURCE, EBCDIC ,N0L 1ST, NOQEU , LOAD, MAP, NO EDIT ,N0ID,N0XREF 
SUBROUTINE FIELD 


ISN 0003 


ISN 000*; 
ISN 0005 

ISN 0006 


ISN 0007 


C 

C 

c 

c 

c 


c . 
c 

c 

c 


CONSTRUCTION OF COEFFICIENTS CMATRIX ELEMENTS') 

IN LINEAR DIFFERENCE EQUATIONS 
SOLUTION BY OVERRELAXATION 

COMMON/BK/IIM,IIP, JJM, J JP , KK . NTOT , I V, JV, II, JJ,M, N, VPC30) , 
1XYZC2080,3).VVC30,20,10),XPC30),XMC1.0),YPC20),YM(10),ZZC10), 
2XX( 40) «YYC30) ,ILX,IUX,KUK,MBC,M8D,VRF,NFPS ,SKPRFL,SKPLST 
CUMMQNVFLD/XC2080,2),COEF(2080,7 ), INDX ( 2 0 8 0 , 6) , SKPCO 
COMMON/CCC/CN,CS.CE,CW,CU,CD,CC,NN,NS,NE.NW,NU,ND 

INTEGER O0,0N/'N0RT'/,0S/'SQUT , /,0E/«EAST*/ 1 OW/*WEST‘/, 

1 OU/'UP '/jOO/'OQWN'/ 

ASSUME ASYMPTOTIC MONOPOLE AT INFINITY 

ALPHAF (UUU)=ABS(UUU/RS) 

NDO=POSITI VE FOR DIAGNOSTIC OUTPUT 


ISN 0008 

ISN 0009 
ISN 0010 
ISN 0011 
ISN 0012 
ISN 0013 
ISN 0014 
ISN 0015 
ISN 0016 
ISN 0017 
ISN 00 1 B 
TSN 0019 
ISN 0021 
ISN 0022 
ISN 0023 

ISN 0024 


N00 = 0 
C 

WRITECM, 1000) 

1000 FORMAT (1H1/18H0FIELD CALCULATION) 

WRITECM, 2000) 

2000 FORMAT C////1X»17HC0EFFICIENT ARRAY) 

X0 = .5*XP Cl V) 

Y0=.5*YP(JV) 

ZGLQ^O. 

DO 600 N=1 , NTOT 

RS=CXYZ(N, 1)-X0)**2 +(XYZ(N,2)-Y0>**2 +XYZ(N,3)**2 
CALL FIND Cl, J.K) 

IFCZZCK).LE.ZOLO.AND.N.GT.l) GO TO 200 
ZOLO=ZZ(K) 

WRITECM, 3000) K.ZZCK) 

3000 FORMAT! // 1 X, 2HZ( , I 2, 2H) = , F6. 3/ 

1 l2X,lH0*l7X,lHSfl7X,lHWfl7X,lHC,17X,lHEtl7X»lHN,17XflHU) 

200 CC = 0. 

C 

C MODIFICATION TO SOLVE HELMHOLTZ EQUATION USING LINEARIZED SPACE 
C CHARGE. HELM = DE B YE-LENGT H-LIKE PARAMETER. CASSUMES POTEN- 
C TIALS ARE DIMENSIONLESS) 

C 


ISN 

0025 


HELM=0 -0 

TSN 

0026 


VOL SQ= 1 . 

ISN 

002 7 


DO 300 MP=1,2 

ISN 

0028 


CALL CNSCMP, C, AREA) 

ISN 

0029 


IF (MP.EQ.l) 00=QN 

TSN 

0031 


IF (MP.EQ.2) 00=0S 

I $N 

0033 


IF (NQO.GT.O) WRITE CM, 888) N , I , J , K , Of) , ARE A, C 

ISN 

0035 

8 B 3 

FORMAT (1X,18HN,I,J,K,00»AREA,C = ,I4,2X,3I3» 1X»A5»1P2E16.4) 

ISN 

0036 


CC=CC+C 

I SN 

0037 


IFCC.GT .0.) GO TO 250 

ISN 

0039 


YYY=XYZCN,2)-Y0 

ISN 

0040 


ALPHA-ALPHAFCYYY) 

ISN 

004 1 


IF CNDO.GT.O ) WRITE (M.999) N , I . J , X , ALP H A 

ISN 

004 3 

999 

FORMAT (1X,14HN,I ,J,K,ALPHA=,I4»2X,313*1PE16.4) 


DATE 



A3i 




ISN 

0044 


CC=CC+AREA*ALPHAFCYYY) 




tsn 

0045 

250 

IFCMP.EQ.l) CN = C 




ISN 

0047 


IFCMP.EQ.2) CS = C 




ISN 

004 9 

300 

CONTINUE 




ISN 

0050 


VOLSQ=VOLSQ*AR£A 




ISN 

005 1 


DO 400 MP=1,2 




1 SN 

0052 


CALL C£W(MP,C, AREA) 




ISN 

00 5 3 


IF (MP.EQ.l) 00 = 0 E 




ISN 

0055 


IF CMP.EQ.2) DO=ON 




ISN 

0057 


IF CNDO.GI.O) WRITE CM, 888) 

N , I , 

J 

, K, 00, ARE A , C 

I SN 

00 59 


CC=CC+C 




I SN 

0060 


IFCC.GT.O.) GO TO 350 




T SN 

0062 


XXX=XYZCN,1)-X0 




ISN 

0063 


ALPHA=ALPHAF(XXX) 




ISN 

0064 


IF CNDO.GI.O) URITE CM.999) 

N.I, 

J 

,K, ALPHA 

ISN 

0066 


CG=CG+AREA*ALPHAFCXXX) 




ISN 

0067 

350 

IFCMP.EQ.l) CE=C 




ISN 

0069 


IFCMP.EQ.2) C W = C 




ISN 

0071 

400 

CONTINUE 




1 SN 

0072 


VOLSQ=VOLSQ*AREA 




ISN 

007 3 


DO 500 M P = 1 , 2 




ISN 

0074 


CALL CUDCMP , C, AREA) 




ISN 

0075 


IF CMP.EQ.l) QO=OU 




ISN 

00 7 7 


IF CMP.EQ.2) 00=00 




ISN 

0079 


IF CNDO.GI.O) WRITE CM, 888) 

N, I , 

J 

,K,00, AREA, C 

ISN 

0 0 61 


CC=CC*C 




ISN 

0082 


IFCC.GT.O..OR.CC.EQ.O..AND. 

MP.EQ 

• 

2) ) GO TO 450 

ISN 

0084 


ALPHA-ALPHAFCXYZCN.3)) 




ISN 

00 85 


IF CNDO.GT.O) WRITE CM.999) 

N.I, 

J 

,K, ALPHA 

ISN 

0087 


CC=CC+AR£A*ALPHAFCXYZCN,3)) 




ISN 

0080 

450 

IFCMP.EQ.l) CU=C 




T SN 

0090 


IFCMP.EQ.2) CD=C 




ISN 

0092 

500 

CONTINUE 




ISN 

0093 


VOL SQ= VOLSQ*AREA 




ISN 

0094 


V Q L = S Q R T CVOLSQ) 




T SN 

0095 


IFCHELH.GT.O.) CC=CC+VQL/HELM**2 



1 SN 

0097 


CALL ARRAY 




ISN 

0090 

600 

r 

CONTINUE 




ISN 

0099 

L 

CALL RELAX 




ISN 

0100 


RETURN 




ISN 

0J 01 


END 






A32 


ISN 

ISN 

ISN 

ISN 


0002 

0003 

0004 

0005 


COMPILER OPTIONS - NAME = MAI 

SOURCE, EBC 
BLOCK DATA 
COMMON/CO/PVOLTS ,X 
REAL PARTI/' IONS* , 
END 


□b/360 FORTRAN H 

N,OPT=02»LINECNT=6Q,SIZE=0000K, 

OIC, NOLIST, NODECK, LOAD, MAP, NOEDIT, NOIOjNOXREF 

MACH,QENST,NN,PARTCL(2>,PART1(2),PART2(2) 

' */,PART2/'ELEC , ,*TR0N*/ 


OATE 



Appendix _B 


Computer Programs: Electric Fields Produced by 

Cloud-to-Ground Lightning Flashes 

The following four pages contain a listing of the computer programs 
written to compute the electric field produced on the ground as a function of 
time and distance from "ground zero" by the charges associated with a cloud- 
to-ground lightning flash. This program was written by Jerry L. Bohannon. 
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T.ITL CLOUD— ro-GROUPD 51MULATII0H 
BATCH 

LAE= STROKE 
0BLG 

IMPLICIT INTEGE R *2 ( I-NJ 
0 1 ME NS i ON RSI<2.1QI- RSIS. <2.10 1 . 

OAT A TP IE/5_ 5 6062E— Ll/^TIP lMy2 ,0E— 7 / 

DATA IC AR0S/'C"y. ITER.K/-M"/ .TI W'r/.IW'N'/ 

DATA I M A/ JT W) L5~-/ . ICY /X".i Cl 6T /. 1 8 E1/X"0 7C 7" /..E8G/ X" 1E10"/ 

OATA IR0/rii011-/.lGR/X^:iC12'/. IYE/X"1D1 3"/. f.6L/X'10l4"/ 

OAT A I A7N/X OE" /• IA7.F/X":0 F" /i.MJLL/ X"0 0' /. IH0M8/X'08~/ 

OATA IBGY/JC IE1 3"/ 

DATA RSIM, 1 >/0.0/.;RS M2 . ll/O.O/.RSIC 2. 10 )/0-0/ 

OATA PIE/ 3- 1415926/ *- 
-OATA RHO/2-OE— 9/ 

CROC T= I. /3- 
URITEC14.1T 
FORMATC'l'T 
00 1000 1=1. 32000 
CCO K=| 

WRITE <14.4> il A7«, iha. I9F.L.IA7F 

FCRMATC 2A2-~ LIGHT KlfflG 30.LT S I MU LA T! I GN-t— .ROI"*2A 2) 

0 WRITEC 14. Ill IBL. TGR 

1 FORMAT < A2." f*E AO *OATA FROM CARDS 'OR T^RHI NAL".-A2 > 

REAOC15.12I : I WHERE 

2 FORMAT < AIT 

3 IFC CWTiERE.EQ. TCAROSMGOTO- 50 

IFC tTIHERE.EO. ITERM) AGCTO. 70 

WRITEC 14. 14* IRQ. IftC 

4 FGRMATC A2."TRY AGAlflT ,A2 > 

GOTO 10 

0 -R EA 0 < 1. 5 l.E ND=9 99 ) T 0. QCL .QSL .VSL ..'VRSJC CRSKUJ). 1= 1. 2T . J=2. 9 » 

1 F0RT<AT<5< F6_ 0 .2 X »/ 8C 2F10 - . C/ T> 

IGOTO 90 

0 WR IT EC 14 . 71 > IMA 

1 FORMAT! A2.' B19TERT FLOAT1 ITG POIRT INITIAL (CONDITIONS F&.0") 
WRITEC14.75 It IBL. ICYT 

5 F0RMAT<A2.' tN IT I'AL HEIGHT KM". A 2 ) 

R£AD< 15.731 tYO 

UR ITEC 14 » 76 ) IBL.ICYt 

6 FORMAT! A2.'CLCUO CHARGE CQUL'.A2>* 

REAO< 15.73# !QCL I 

'WRITE < 14 . 72 ) I8L. ICYT 

2 FCRMATC A2 ."STEPPED LEADS R CHARGE CCUL'.A2) 

REAOC 15.7 3* :QSL 

3 FORMAT < F6.0 I 

WRIT E< 14. 74T IBL.ICYT 

i- FQRMAT<A2."STEPPE0 LEADER V8LCCITY E5 M/S".A2> 

R EAO (15.731 IVSL t 
‘WRITEC 14.77):. IBL. ICT 

7 FCRTIAK A2 •" RETURN STROKE! VELOCITY E7 M/S".A2) 

READ! 15. 73* :VRS ‘ 

I URITE<14.aO) IBL.ICYl 

) FCRMATC A2 ENTER 3 -TICES CMS) AND CURRENTS (JCAMF) TO DEFINE THE”!/ 
«" RETURN STROKE 2F10 .C'/t' i'.!9>. " ! A2 ) 

00 32 J=2.9 

REAOU5.31) RSI <7. J>. RSI < 2- J) 

IFCRSIC 1. J». LT. 3. » IGCTO. 73 

CCNTINUE B2 



31 :FORMAT*2F10.. 0 ) 

90 VSLS=VSL 

RSI Cl. 1 0 )=RS I *1. 9 > 

00 10 02 J=l. 10 
RSIS* l.J)=RS I Cl. J* 

RSI5*2.-J)=RS 1*2. J» 

RSICl. JI=RSl C li -11/1000. 

RSK2. J l = RSI < 2m— 1 ) » 1 0 0 0 . 

L0LO2 CONTINUE 
QCL=-GCL 
QSL=-OSL 

sycs=iyo 

RAO=< 0. 75 *A. 8SCQCL I /"R1H0/ Pi I El+^CRQOT 
00 1005 1=2.9 
A=RSI C 2 • I ) 

B-RSI <2. I— 1 1 
C=RSI< 2.1+1! 

rF<A^GT.B.ANO.A.:GT.C 1 IOT = I 
l C 0 5 ICCNTINUE 
VRSS=VRS 

LOO VSL=— VSL * 1- 0 £5 

YO=TO*1 00 C. 3 
>VRS=VRS*1.CE 7 
OISL=I.0E-4 * 

105 WR I T EC 14 • 11 0 ! IBL.ICT 

L1C FORMAT* A2."W«AT IS RA0UfS'.A2» 

■READ* 15.7 3* O 
URITEC13. Ill I IRD.lffGY 
111 FQR-MAT* 2 A2/"-L"> 

00 1 00 1 1=1. 320 OC 
.001 K=1 

WRITE* 1 3 * 149 » 

1 49 FORMAT* 1X."S:I UNITS"//! 

.50 'WRt T£ * 1 3. 151<! ; YCLS.QCL.QSL »VSL5» VRSS.RSI S'. 0 

L5 1 FORMAT* IX." «E IGHT= **. F7_;0." KM" / 1 X . "<J— CL OUD = ". F7.1." C" / 

SlX."0-LEA0ERf= ".F6-1." C.' /I X . "V- LE AOER= ".F6-1." E5 M/tS"/lX. 
S "V-RETURN= "•» F 6. 1. " E7 M/S' /l X." RETURN STROKE MS. KAMP'/ 

Si :<2F10.4/»/‘y lX."RAaitiS= ".F6.0." Ml"///)' 

WRITE*! 3. 152 » IRG. IRD 
52 FORMAT < IX. 

SA2.aX."T'.15X."E".15X."Q;',l&X."H". A2» 

T = 0_ C 

SLR Y 2= 1 - 0/t Of* D* YQ+YHJ 
SLRY=SQRT*SL RY2 » 

• YC=Y.0+R AO 

SLRQC2=1. 0/CD+D+ YC*YC> 

SLRQCL=SQRT *ISLR QC2F+.SLR0 C2 
01 = 1 . 0/0 
X = YQ 

EMAX-0.9 
G O CONTI l NUE 

SLRX2=1 ,.0/C a*D+x+x» 

SLRX=SQRT*SI_RX2 > 

SLR X 3 2= SL RX +‘SLR X2 

E=QSL/TPIE/rO**SLRX-SLRr.)-*SL'R CCL* YC/TPIE* * QCL-“QSL+< 1.-X/YOI » 
IF* ABS* E ) -GT . ABS* EMAX ) ) ,EMAX=E 
IF* ABS* E ) .LTf.5.0 E4) GOTO 211 
WR1TE<13.210 »T.E.X 
LG FORMAT* F16.7. FI 6.1. 16X. F 16 - 1 ) 


B3 



215 FORJIATt F16-7. F 1*6.0, F16..5 ,F 16 .1 . I I C.F 16 . 7, F 16 . ? ) 
211 IFTABSt E> .LTT.5. C£4. OR-X-.GT.O .SE3 > DTSL=1.0E-3 
T=T-+DTSL 

I Ft X- LT . 50- > GOTO 50) 

X=YO+VSL *T 

IFTX.LT. J.O I GOTO 50 ) 

DTSL=1 . 0£— 4 
GOTO 200 
30 0 iCCNTINUE 
T=T-OTSL 
UR ITE<13»501 I 
301 F CR MAT 1 1 X*" "I 
ESL-E 

QRC=QSL-QCL I 
PL=-QSL/YO 
SLRY.03= VC *SL ROC Li 
■KRNT= 1 
R 1=0-0 
i T R= 0 - 0 
KGtJ>=0 
0 = 0-0 

ilc CONTINUE 

CALL; CURE NT ttR SI ,0,0*1, TR,IR J.KiRAT.XOLD* 

I F< R I . L £ . L- 0 > GOTO 6:. 

T=Tn»-OT 

Y=VR3*TR 

IFtY.GT.YGI (GOTO 522 

p=oy r 

SLRYR=1 -C /SaRTCa*0*Y.* Y» 

E=ESL-*P»><DI-SLRYR J/T-PIE 
I F< ABSt E > -GT*. ABSt EMOJX I » ,EMX=E 
IFtTQT.LT.KRrNT.ANO. ABSTE>.lT.5.CEO GOTO* 51" 
WRITE* 13.215 ) T„E,CU.Y.KR.NT ,)TR.RI 
GOTO 510 

22 'WRITE! 1 3. 50 1 ) 

20 CONTINUE 
P=Q/'YG 

IFtP.GT.PL* !GOT[r 572 
E=ESL-*P* < 01— *SLR Y J/TPIE 
IF< A8S t E ) .GT . ABSt EWA r X )> .EM X=E 
IF< 1QT.LT.XRNT. ANO.ABStE I.LT.5.0E4) GOTO 521 
WRITE (I 3,215 I T»’£«0,(YG,X!RM,!TR» R I 
2 1 CBLV. CURENTtIRSI ,-Q,OT.TR,lR I.XIRAT.XOLD* 

IFtRI.LE .0.0 ) GOTO 6.0 

T=T-*DT 

GOTO 520 

f 2 iWRITE C 1 3, 501 I 
TO ICGNTINUE 
QRS=Q>QSL 

IFtGRS.GT.QRfO GOTO; 600 

E=£SL-*PL* <0 T-SLRlY 1/TPIE*)QRS* SLR Y 03 / TP [E i 
I Ft A8S< E > ..GT -A8S4 EMA'X ) J .E MAJ6= E 
IF< IOT-LT.KJTNT. ANO. ABS< E) XT.5.0E4 » GOTO 571 
•WRITE 1 13.215 > T..E.Q.IYC, K!RM,iTR. RI 
1 CH.il CUR£NTtlRSI«0,OT.TR,lR I.K'RAT . KCLOJ 
IFtR I.LE.C.O > GOTO 6. 0 
T=T^OT 
GOTO 570 

j WRI TE < l 3 , 59 9 I Q ,EMAX 
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El 2- 4 


V/MT » 


99 FORMAT!// 1X„ "CR T= ',F10.,4 f * C".5X.'EHAX= 

‘WRI TE ! 14 » 60 1 I 1 BG, IBL,IGR»IKEL 

1 FORMAT! 2A2." DO TOU -VANT IAAOT HER R A D [US' . 2 A2 > 

2 READ ! 15,121 ilAO 
IF<IAD.EO.ir > GOTO 1.5 
IF! IA0-EQ.I VI GOTO 650 

'WRITE! 14. 141 IRD.IGR; 

GOTO 602 

50 UR1 TE ! 14, 65 1- ) 1 8t_ . I GTJ 

51 FORMAT! A2 ,“ 0-0 YOU WANT AINCTHER EVENT*, A2) 

52 REAO!15,12» IE 

IF! tE-EO. ir >: GOTO 13 
IF! IE. EO- IN* GO tO 95.J 
WRITE! 14, 141 IRD.IGR; 

GOTO 65 2 

99 WRITE! 14, 99 8 > IRD 

98 FORMAT! A2,~K0 MORE CAROS' I 

50 WRI TE! 1 4* 95 lrf I fTA, ICR 

51 FORMAT! A2.~ END OF PRUGR A M", A 2 > 

STOP 1 

ENO 

L A8= CURENT 
BBUG 

Subroot ine curent<rsi,q,;ot, t;r,ri ,kirkt,kold» 

IMPLICIT INTEGER -2 ! I-M* 

DIMENSION RS 112,1') . 

IFIKOLD. EQ-KRNTIi GOTO 50 
TA0=RSI!1.KR!NT«- LF-RS 1 < 1 ,IKRNT ) 

IFTTAU-LE.O-C) GOTO -ICO 
IF!TAl).LE-l- OE-5 I DT=0.5‘E-6 

IF!TAU.GT.U CE-5- AMa.TAIC.tE. l-OE-4 » 0T= 1 . E-5 
IF! TAD. GT .1,0 E— 4 - ANCI- 7 Alt . L£ . 1 ,0 E— 3 » 0T=1. £-4 
IFtTAU.GT.l- 0E-3-AND.TALK.LE. l.OE-2 ) DT=1. E-3 
IF!TAU.GT.l- OE-2.ANO.TAir.LE .’1 .OE-1 » DT=l.-E-2 
IFXTAU. GT.l. CE-l 1 DT.= C.0-25 
OEL 1= IRS I !2-» KRN T?* IF— RS l <'2 WCRNT) F/TAU 
IF!G.EQ.O.O» R-1 1 = 3 . 0 
TRR=TR+DT 

IF!TRR.GT.RSUl,KRNTi> U *. CT=RSI i 1, KRNT*- 1F-TR 
RI=Tri>OELI*OT 
!RI2=RI 

Q=0-*DT * tR I2-«- R III/ 2. 

TR=TR-*DT 

■kcld=k:rnt 

IF! IR. GE . RSI ! 1, KR NT-*- 1 M )KRNT=KRNT+ 1 

R 1 1= R I 2 

RETURN 

lu CCN71INUE 
R 1=0-0 

return 

ENO 

;£f»D 
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Appendix C_ 

Computer Output Listing: Cloud-to-Ground Lightning Flash Density 

The following seven pages are the computed output from the program that 
calculates the lightning flash density (cloud-to-ground) from the monthly 
thunderstorm days using the Pierce Conversion. This program, written by Jerry 
L. Bohannon, uses the Normals, Means and Extremes data from "Local Climatolog- 
ical Data — Annual Summaries for 1977" published by the National Oceanic and 
Atmospheric Administration, Environmental Data Service, Asheville, North 
Carolina (available also on magnetic tape). 


Cl 


Cl] 

123 

[3] 

[43 

C53 

[£] 

m 

b] 

i?3 
[10] 
ill] 
[ 12 ] 
[13] 
*■143 
C153 
C16j 
[17 j 
[133 
[1?] 
[203 
[213 
[223 
[233 
[243 
[25] 
[261 
[27] 
[283 
[ 2 93 
[301 
[31] 
C 321 
[33 ] 
[343 

lo OJ 

[363 

[373 

[33] 

[3?j 

[403 

[41] 

[423 

[433 


7C3MFUF0T 

7 AARj-WID CO-ifUFOT Hi 

ft THIS r UNCTION COMPUTES THE ELECTRIC POTENTIAL IN 3 RESIGN A R ROUND 
ft ONE BILLBOARD GF THE RECTENNA, THE MEASUREMENT AREA STARTS METERS 

| FROM THE LEFT HAHD EDGE OF THE REGIE HA AND EXTENDS TO THE RIGHT 1 SC 1 MET 
j THE BOTTOM OF THE MEASUREMENT AREA IS AT GROUND LEVEL, SMILE THE TO? 
a IS 'UP* METERS HIGH, 

H THE RESOLUTION IS CONTROLLED BY THE ARGUMENTS OF THE FUNCTION, THE 
ft LEFT ARGUMENT SPECIFIES THE NUMBER OF COLUMNS IN THE OUTPUT, THE RIGHT 
ft ARGUMENT IS THE HUMBER OF ROWS , 

il THE FORMAT OF THE OUTPUT IS AN ARRAY OF NUMBERS IN SCIENTIFIC NOTATION 
ft WITH ONLY ONE SIGNIFICANT DIGIT PRINTED, 

?OTi-(HI,SID)fQ 

CF-F'H 

Gi-Rf-1 

LOCPHJH1p(Q-J)l‘JPfHI-J 
Rf-1 • 

LOOPL *L1f J1 , %f( R-i)XACf«D-l 
P0T[0;R](.L) FIELD H 1 
RfR+1 

4(^1 WID) /LOQPL 
GfOj-1 

4(Q<HI)/L00FH 

TSYf-TSY+l 

M f-lt3?u TS ■ 

■THIS IS RUN HUMBER ' , ( t TRY) .DATE 

'THE CALCULATED VALUES CF THE ELEOTRIC POTENT I A*. , IN vCLTS, ARE SHOWN BELOW 


AARfFOTj-gpOT 

5/a«vr 2 0l] ' 
aFFpio 

1 “HE VECTOR OF lIHE CHARGES USED I5 t## CC‘ULGr-c‘3 f SR HE^ER, 1 

‘THE SUM OF 7FE LIME CHARGES 15 1 f (f+/LA) t 1 COULOMBS PER METER,' 

1 THE TOP OF THE MEASUREMENT ARRAY IS METERS HIGH, 1 

‘THE RIGHT EDGE OF THE ARRAY 15 ' METERS F*C.M THE FIRST DILLSCflRI- , 1 

•THERE ARE 1 f (f j , 1 COLUMNS PER: METER, A’-fS ; j ■, fUF-T ; » ? TORS PER .i.STR 

1 5:UH MG, 1 ,( f TRY) , DATS 

QFF>1 

w/3 fiV £201J 

4(SI3HrQ)/0 

1 1 n — Artp'Af 3z.LGW iHuWS THt SIuN Or i.ACH OF THE i tLMBZRS IN THE ABOVE ARRAY 

i i 


XFOT 

'THIS IS RUH HUMBER ' , (fTRY) ,DATE 
7 
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"PROTECT [OjV 
7 MID PROTECT HI 

[1] „ THIS FUHCTIOH COMPUTES THE ELECTRIC POTENTIAL VA A REGION OF SPACE 

[2] s DUE TO A CHARGED MIRE LOCATED SOME FIXED PERPENDICULAR DISTANCE -ECM 

[3] D THE TtlP OF EACH BILLBOARD OF THE RECTEHHA, THE MEASUREMENT AREA 15 

[43 <| EXACTLY THE SAME A5 THAT USED IN { <C3HFUP3T> > , A5 WITH ( <COj*FUPOT> > 

[53 # THE RESOLUTION IS DETERMINED BY THE ARGUMENTS OF THE FUHCTIOH, 

[6] ft the FUHCTIOH DOES HOT PSIHT ANY OUTPUT, THE OUTPUT IS CONTAINED II! 
[7jj 0 THE VARIABLE, (<PROT>), THIS VARIABLE MILL HAVE THE SAME DIMEHTI0N3 A 

[8], ft {<POT>>, THE VARIABLE CONTAINING THE OUTPUT FROM < (CSMPUPOT) » , 

[?] PROTp(HI,MID)fO 
CIO] OlfftlH 

[113 A { ; LOOPH) } COMPUTES ALL OF THE VERTICAL INDICES, 

[12] looph;h2p(Q1-1)xUPshi-i 
C13] Rl«-1 

[14] ft <<LCOPL>) COMPUTES THE HORIZONTAL INDICES AND CALLS <<F!ELB»)!, 

[153 L00PL*L2(-31.?it(Rl-l)xflCi.5(iD-i 

[16] PRCT[si;si3pL2 field* «2 

[17] sipsin 

[13] H(RUi«D)/LOOFL 

[19] ait-ai+i 

[20] H(31> HI) /LOOPH 

[213 PROTpePROT 

[22] TRYIpTSYI+1 

[233 'THIS IS RUN NUMBER 1 , (f TRY'D j 1 OP PROTECT DATE 

[24] DTP HO 

[25] 'THE PROTECTING MIR: E IS LOCATED 1 , lt"“'l) I ' METERS FRONT THE' 

[26] 'LEFT EDGE OF THE ASSY, AND 1 , ( fSl jX-Jl ; , : METERS FROM THE BOTTOM,' 

7 


C3 


7FIEL0 [03? 

7 UHLi-L FIELD H 

S THIS FUNCTION CONFUTES THE ELECTRIC POTEHTIAL AT s;j f rOIHTj £LjH) 
[2] fl IH.THE SPACE ARCUMD THE A SR A f OF FIVE 3ILL33ASDS, 

cn i-iFiS.yjx'i+if > 

[43 

[53 s«-3e32r9 

[33 Af-L-XJ 
[7] U-t 
[33 uif-H/O 

[9] BSLOOP;HM{.((H-Sx:<J)jt2)T«A}.(A-LI[I])t2 
CIO] DMf((H+SxXJ)*2)+ Hft 
[113 UI[I]F+/-(LflP32xEO)X»()lh>SM)j0,5 
E123 IH+1 
[133 -»( I<H+l)/SBl_OOF- 
[143 uhlf(-100003xH)+t/«i 
? 


7FIELDW [Q]7 
7 PfL FIELDB H 


[13 

ft THIS FUHCTIGH COMPUTES THE ELECTRIC POTEHTIAL AT ANY POINT, (L,H), 

[23 

ft DUE TO THE CHARGED PROTECTION WIRE ABOVE 

THE BILLBOARD, THIS WISE IS 

[33 

ft ASSUMED TO BE PARALLEL TO THE BILLBOARD 

AHD LOCATED A PERPENDICULAR 

[43 

ft DISTANCE, ((SPACE)), FROM THE TOP OF THE 

BILLBOARD, 

[33 

ft THE CHARGE OH THE WIRE IS ;('_*>), 


[33 

'-ilf-15, ?3*'1+|H 


[73 

Mf, L Wf- , LW 


[33 

LOHGPl2,24r23THTfip-3.3(Pl2,24)X5PACE 


C93 

:<j1p,;;jj>lciigx20Thta+02p9 


[103 

s lP33(02f-9)pTHTA 


[113 

AipL-;;ji 


[123 

ip 1 


[13.3 

uitfH^g 


[143 

loop ;-::MiD{(H-5ixxj2)*2)+ Ha :P(Al-Lii [i]) »2 


[153 

0«lt-( C”+51 xHJ1)*2J+n«1 


[133 

“1 1 L : ](■ +/- ii-il rvZ'i Z J)- 1 ( N'H ( f I M i ; ,g • 5 


[173 

If-Iti 


[133 

+(I(H*l)/LOOP 


[1?3 

ppp/ui! 
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5T A TZ 

STATION 

THUNJc^ST C R M 

GA C ‘JND 'S7-7 IK. 



D A Y j 

DENS I T Y 



(NO »/Y?Art ) 

( nC./Yp. /KM 

AL 

3 I 9 Vi I N Grt A M 

5 3 ■> 7 \ 

.13,37 

At_ 

HUNTS V Il_L £ 

:v. 2 0 

: 3 , 3 a 

A L. 

M3 3 I LZ 

79,75 

27.07 

AL 

MONTGOMERY 

6 2,12 

1 ^3 * 4- 


AK 

AN CHOP AGE 

1.1-2 

0.39 

AK 

ANNETTE 

1 , 43 

0.60 

A K 

BARRC 

U. 07 

0.3ft 

A K 

BARTER ISLAND 

0.26 

0.17 

AK 

BE 7 HE L 

1 , p9 

0.5 1 

A K 

BETTLE5 

4.67 

0.73 

AK 

BIG DELTA 

2,5 0 

0.41 

AK 

COLD BAY 

0. OS 

0.04 

AK 

r A I Pu ANK S 

3.12 

0.33 

AK 

GUlK ANA 

4, 7 0 

0.76 

A K 

HOME 9 

0.33 

0.15 

AK 

JUNEAU 

0, 32 

0.24 

A K 

KING SALMON 

1,31 

0 . 48 

A K 

KOOI AK 

0, 29 

0.09 

AK 

KQT ZLBUE 

0, o2 

0.24 

A K 

MC G 9 A TH 

6. 77 

1.05 

AK 

NOME 

0.42 

0.19 

AK 

ST PAUL ISLAND 

0. 3 6 

0.05 

A K 

SHEYMA I 3LANJ 

0,15 

0.10 

A K 

SUMM I T 

5,0 0 

0,32 

A K 

TALKEETN A 

4,30 

0.77 

A K 

UNAK1_EcT 

l , ?0 

0.46 

AK 

YAKUT AT 

1 . 59 

0.60 


AZ 

FLAGSTAFF 

50. 53 

15.37 

A Z 

PHOEN IX 

2 3.03 

4.60 

A Z 

TUCSON 

3 9.34 

13.29 

A Z 

w I N5L C W 

36, 34 

9.3 3 

A Z 

YUMA 

7.26 

1.4 3 


AR 

FC P T SMITH 

o7 .06 

11.64 

A PK 

LI TTLE ROCK 

5 6.97 

11.51 

A S 

PAGO PAGO 

26. 09 

3.73 


C A 

OAK EP S F I - 1_ D 

2,3 0 

0,9^ 

C A 

6 13 HQ P 

1 3, 24 

2,21 

C A 

BLUE CANYON 

11.33 

1 • 7 

CA 

eUREK A 

4,56 

1.23 

C A 

FRESNO 

5, 43 

1.37 

C A 

LONG BEACH 

3. 71 

1.10 

C A 

LOS ANGLES (CITY) 

6,21 

1.43 

C A 

LOS ANGLES (L AX ) 

O.ol 

1.10 

C A 

MOUNT SHASTA 

1 3, 2 7 

2.16 

CA 

OAKLAND 

2,23 

0.37 

C A 

RED dL UFF 

9. ^0 

1 • 32 

C A 

SACRAMENTO 

4,76 

1 • 27 
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I I II II 


5T ATE 

5T AT I ON 

T HUNDE- STORM 
DAYS 

( NO. /YEAR ) 

G- 0 UNO -> S T-. I K. 
D-NS ITY 
( N3./YS . /AM 

C A 

SANG3EP0 

a, 22 

1,14 

CA 

SAM 0 I =.GJ 

cl » 7 0 

0.56 

C A 

SAN FRANCISCO (CITY) 

2,05 

0.36 

C A 

5AM FRANu I SCO ( 3F 0 ) 

2, 1 2 

0 1 

CA 

S T Q CK T ON 

3 . l 1 

l . 0 1 

CA 

SANTA MARIA 

2 . 32 

G • 6 9 

cc 

ALAMOSA 

44.A2 

12.92 

cc 

COLORADO SPPIN-JS 

5 9 ■ 6 9 

22.43 

CO 

DENVER 

4 1,33 

11.02 

CO 

GRAND JUNCTION 

34, 32 

6.38 

CO 

P‘jEfct-0 

a 0 * 3 2 

10.59 

C T 

8RI0GEP0ST 

■21.57 

3.50 

CT 

HARTFORD 

'22, 30 

3.62 

OE 

WILMINGTON 

31.03 

5.7 3 


DC 

WASHINGTON 

(OCA ) 

29, 07 

5.13 

DC 

WASH! NGT ON 

( I AD ) 

27.13 

4- • 6 0 


fl 

APALACH I COL A 

70.19 

2 2 . 9 9 

fl 

DAYTONA BEACH 

7 9.61 

29 . 25 

fl 

FORT MYERS 

9 4,67 

4 7,04 

fl 

JACK SONV ILLS 

o 3 . 9 4 

20.2c 

F L 

KEY WEST 

6 2.66 

15.38 

FL 

lakeland 

99. 5 0 

4 3.56 

FL 

MIA K I 

74,04 

26. 37 

fl 

ORLANOO 

8 1.21 

32.79 

FL 

ORLANDO (MC COY AF3) 

^3,62 

30. 37 

FL 

PENSACOLA 

74,13 

22,90 

FL 

TAlLAHAS see 

56,3 7 

35 , 05 

FL 

TAMPA 

83. 19 

40.90 

FL 

WEST PAlM BEACH 

7S.63 

23.69 


G A 

A T m2 N S 

5 1.52 

1 3 .00 

G A 

ATLAN TA 

50 a 1 9 

11.57 

G A 

AUGUSTA 

5 5.15 

15.41 

G A 

COLUMBUS 

56.7 1 

15.61 

GA 

MACON 

56, 33 

la.43 

G A 

ROME 

6 1.42 

15.37 

G A 

SAVANNAH 

64,33 

20,62 


G U 

TAGUA C 

27.0 5 

4.79 

H I 

HILO 

3. 75 

1.67 

H I 

HONOLULU 

7, 07 

1,43 

H I 

K.AHULU I 

4. 95 

1.20 

H I 

LI HUI 

3. 5 1 

1 .64 
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ST AT = 

STATION 

thunder storm 

jri 5UNG~»ST A> I 



DAYS 

3 z.N S I T Y 



{NC./YEAR ) 

( iN - a / Y P »/K. 

I C 

83 I 3E 

I A , 

2a 34 

I D 

L~ w IS TON 

I 3. 73 

2.45 

I 0 

POCATELLO 

03.11 

4 .02 

I L 

CAR 1 3 

52.77 

10.9a 

I L 

CHICAGO ( MID'* AH' J 

A 0 . 5 a 

7.53 

I L 

CHICAGO { O* HARE ) 

3 5. 4-2 

O ^ 7 d 

I L. 

MOL IN £ 

47. 36 

10.01 

I L 

•^ECRI A 

Ay i 44 

10.25 

IL 

ROCKFORD 

a 2 • 19 

5.33 

I l_ 

5PR IN GF I EL C 

5C > JO 

10.73 


I N 

SV AN 3 V ILLS 

45.73 

3.87 

I N 

FORT WAYNE 

4 1.00 

7,37 

I N 

INO llANAPQLl S 

4 4-. 6 9 

5.37 

I N 

soutU SEND 

4 2.39 

3.64 

I A 

BURL I N GTGN 

50,55 

11.06 

I A 

GE S MOINES 

4 9.73 

11.22 

I A 

QU8UGUE 

4 4. '9 o 

9 . 29 

I A 

SIOUX city 

45. 38 

10.49 

I A 

WATERLOO 

4 1.70 

8 . a 1 


K. S 

CONGO 3D I A 

5 3.43 

15.71 

K. S 

DODGE CITY 

5 3 , 03 

14.33 

Is. S 

GOGOL AND 

4 3.74 

1 3 . o 9 

R S 

TOPEKA 

6 7.58 

14.14 

K S 

nICHITA 

55, 29 

13.25 


< Y 

LEXINGTON 

46,7 6 

10.22 

K Y 

LGU 13 V ILLE 

4 5,40 

9.13 


L A 

ALEXANDH I A 

53.07 

16*96 

i_ A 

BATON ROUGE 

7 0.46 

20.07 

L A 

lake Charu.ES 

7 6.38 

22 . 59 

L A 

NE * ORLEANS 

6 6.43 

20 . 38 

L A 

SHREVEPORT 

5 4,15 

10.21 

ME 

CASI30U 

20.33 

3.67 

M £ 

PQF.Tl and 

l 8. 05 

2 . 9 6 

M D 

bALT I MCRE 

23,44 

5.10 


M A 

30STGN 

1 9, 33 

3.14 

M A 

NANTU CKET 

20. 27 

3.09 

M A 

wQRCc. 5 Tc.R 

2 1,27 

3.51 
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5T A TE 

s'a n jin 

T tiuNOE Tj, O R,-t 
DAYS 

(NU,/YEAR ) 

_/U ND 5 7 i- I 
DENS! i Y 
( NC ,/ YR . /.< 

M I 

AL P £M A 

33. 29 

6.25 

M I 

D5TSCIT (OTT) 

32>02 

5.67 

M I 

OST-iOtT (OTA) 

3 3, 2 = 

D • o 7 

M I 

FLINT 

3 3, 03 

5.97 

M I 

GRAND SAP IDS 

36.7 1 

O . a 

M 1 

MdJGn TON LAKE 

3 do da 

7.6 5 

M I 

LA NS I N o 

04,17 

6.18 

M I 

MAR DUE TTE 

28.07 

5.0 3 

M I 

MUSKEGON 

37, ia 

6,83 

M I 

SAUlT STS MARIE 

29. 4a- 

5.22 


MN 

CJLUTH 

3 4,66 

7, 3a 

M N 

INTERNA i ICNAL rALLS 

3 1.42 

6.67 

M N 

MI NNE APOL I S 

36.79 

7.4 1 

MN 

ROCrlc 3 TE R 

4 1.00 

6,3 2 

MN 

SAINT CLOUD 

3 5.76 

7 , d4 


M S 

JACKSCN 

6 5.79 

16.30 

MS 

MLR I DAN 

5 6.59 

1 3 . 3 1 


MC 

CQLUMS I A 

51,50 

10.40 

M u 

KANSAS C 1TY ( MCI ) 

5 1 . 20 

11.39 

MG 

KANSAS CITY (MKC) 

49,66 

10.59 

M 0 

SAINT jCSEPh 

5 6.35 

13.76 

MC 

ST. LOUIS 

4 4.55 

. 6.61 

M G 

3PR INGFIELO 

58.00 

13.00 


M T 

81 LlI NSS 

2 8, 79 

6.00 

MT 

GL A5G0 w 

27, 1 1 

5.30 

M T 

GREAT FALLS 

2 o » 60 

5.17 

M T 

HA VRE 

21,60 

3,86 

MT 

HELEN A 

33.6 1 

5,32 

MT 

KALI SPELL 

22,75 

3,5 0 

M T 

MILES CITY 

26, 43 

c « 0 2 

MT 

MISSCULA 

23, 61 

4.3b 


N E 

GRAND ISLAND 

4 7 . 9 9 

11.76 

N E 

LINCOLN (APT) 

46,33 

1 0 . 77 

N E 

LINCOLN (CITY ) 

4 9 , 3 a 

11.99 

N E 

NO R FO L K 

50 ,20 

13.11 

N c 

NORTH PLATTE 

45,92 

11, 55 

NE 

OMAHA (CITY) 

4 0, 50 

3. CO 

N d 

OMAHA ( EPPLY FIELD) 

4 3.60 

l 1 , 2 f 5 

N c. 

SCOTT 56L UF F 

4 5, 56 

11,92 

ME 

VALENT InE 

45, 22 

12.75 


N V 

ELK 0 

20, 72 

3.47 

N V 

ELY 

3 2 « 0 0 

c • 7 o 

N V 

LAS VEGaS 

1 A, *7 

2 • a 5 

NV 

RENO 

13, 5 & 

2.06 

N V 

Wl NNE MUCCA 

l ^ » 3 o 

2,24 
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0 


5T ATE 

3T AT I ON 

ThuNO--? 5 T j R .*• 

3-' :un 0->5 r?> 



OATS 

Oc'i 5 I TY 



( NJ. /YEA A ) 

( NCj/YA.A 

N H 

CONCORD 

2 J , 4 7 

3 a 4 9 

M H 

M T A 3ni NGTC-N 

1 =>, 33 

2,7a 


N J 

ATLANTIC CITY 

23, AT 

A , 3c 

N J 

N£ ,* ARK 

25. A 7 

A , A G 

NJ 

TRENTON 

3 3,22 

0,35 


NM 

ALBOQU EriOU E 

42,34 

11.13 

N M 

ClA YTON 

54.11 

17.33 

NM 

RG SiicLL 

3 2 , o 0 

6.30 


N Y 

ALBAN Y 

27, c 4 

5.20 

NY 

81 NGHAMT3N 

3 1. a 2 

5.94 

N Y 

6UF r AL <3 

30. 74 

5.13 

N Y 

NE M YORK (Cl TY ) 

19.A7 

3.16 

NY 

NE* YORK (JFK) 

22, 32 

3. 56 

N Y 

N = N YORK (LA GUAhOIA) 

24,24 

4.31 

NY 

RQ CrlcSTiR 

29. 2a 

5.21 

NY 

S Y RACU 5E 

29.3? 

5.43 

NC 

A 3 H E V I LL E 
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